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THE TWO-DIMENSIONAL IRROTATIONAL FLOW OF 
A COMPRESSIBLE FLUID IN THE ACUTE REGION 
MADE BY TWO RECTILINEAR WALLS 


By J. WILLIAMS (£zeter) 
[Received 2 July 1948] 


1. Summary 

An exact solution is derived for the two-dimensional irrotational 
flow of a non-viscous compressible fluid in the acute region contained 
by two inclined rectilinear walls. In the case when these walls are 
mutually at right angles, tables are given to illustrate the deviations 
in fluid velocities from the classical incompressible case. 


2. Introduction 

I shall assume that the pressure-density relation follows the 
adiabatic law. Also, following Temple and Yarwood (1), I introduce 
a non-dimensional speed variable in +, where + = q*/¢?, and where q 
denotes the fluid speed at any point and g the theoretical maximum 
speed for the motion. Hence g? = 28a, where B = 1/(y—1), and ay 
is the acoustic speed in the locality of a point of fluid stagnation. 

The solutions of the hodograph equations of motion have already 
been given by Chaplygin (2), a typical solution being 


(7, 0) = R,,,(z/7,)'"F,,(7)sin m(9+-€,,) 
6 km 
(7, 0) = —*Be rye (=) F(x) }eos m(0-+ey) 


m T\\T, 


i(m—1) 
7 qz(r, 6) = Rk, —t) (=) x 


Ty 


< (F'n (r)-+-mF,, (7) é —iw al TF,(7) ew 


i l—m 1+m 





provided that m + 0, +1, and where 
w = (m—1)0+me,,, F_(r) = F(A, B;m+1;7), 
A+B=m-—B, 2AB = —m(m-+1)B. 


Again, o, ¢, 2 = x+iy, 6 have their usual meanings, R,, and e,, are 


~ 


[Quart. Journ. of Math. (Oxford), Vol. 20, Sept. 1949, pp. 129-34] 
3695.20 K 
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constants, and 7, is some standard value of r. A prime attached to 
F,, denotes differentiation in r. 

The ‘corresponding solution’ in the case when the fluid is assumed 
incompressible can be deduced from (1) by making g > 00 in such a 


way that r'¢ > q and 73g > U (say), givin 
5 q>q 19 y), giving 


$(q,0) = R,,(q/U)™ sin m(O+-€,») 
$(q,9) = —R,,(q/U)™ cos m(6+-e,,) . = 
U2(q, 0) = R,,(q/U)"™-“K{m/(1—m)}e-”  (m #0, +1) 


It is quite straightforward to show that the equation of the 
envelope of the curves of constant — corresponding to (1) ist 


(1—r)(27F',, +mF,,)? tan*{m(8-+-e,,)} = {(28+1)r—l}(mF,)?, (3) 
and, consequently, lies entirely in the transonic and supersonic 
regions of flow. Furthermore, in the manner of Temple and Yar- 
wood (1), it can be shown that on the envelope (which can be 
regarded as the boundary separating regions of the physical plane 
in which 7, 6, ¢, and % are either single- or multi-valued functions 
of x and y) the streamline possesses a singular point, which takes 
the form of a cusp. Theoretically, in the case of those streamlines 
which have no singularity in crossing the envelope, the fluid speed 
can attain both subsonic and supersonic magnitudes without contra- 
vening the conditions of continuous flow. Consequently, the highest 
velocity attainable without the formation of discontinuities, i.e. 
shock-waves, is dictated by that streamline which possesses a cusp 
at the lowest value of 7 in excess of 1/(28+-1), the critical value. 

In the case of the solution represented by (1), which is not treated 
in (1), I quote the final equations which determine the values of 7 
corresponding to the cusps as either 

(i) mF,,+27F), 
or (ii) mBF 2 ml 28+-1)r—1]+-(28+-1)7}+ 
+2F,, F’,{2m(1—z7)[ (28+ 1)r—1]+-B(28-+- 1)7?} 
+4F)?7(1—7)[(28+1)r—1] = 0. 

These equations were obtained by expressing the condition that the 


equation for 4? found by eliminating @ between (1) and (3) should 
possess double roots. 


(4) 


+ This result is not quoted in (1). 
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3. The two-dimensional flow in the acute region made by 
two rectilinear walls 
Consider the complex potential 
Q = ¢+if = —U" (n> 2), (5) 
which represents the two-dimensional irrotational flow of an incom- 
pressible fluid between two walls as rigid boundaries and containing 
an acute angle « = 7/n radians. We can express z, and hence 7 and 
¢, as functions of g and @ as follows: 
dQ 
dz 
so that z= (q/nU)Hm—Detr-Oin—1) 


= ge- — —nUz"-1, 


since O<argz<ac}hn, r>0Da. (6) 


Consequently, in the region 0 < argz < a of the physical plane, z is 
a one-valued function of g and 6. Referring to (2), we find that this 
solution corresponds to 

m = n/(n—1), R 
from which we write 


h(q, 8) = —U(q/nU)™— sin[ (nd —z)/(n—1)] 


= —U(n)-™n-9, Em = —7/N, 


m 


7 
$(q,0) = U(g/nU)"*"-® cos{ (nd—z)|(n—1)] ” 


The solutions ‘corresponding’ to these in the case of compressible 
flow are written down from (1) as 
—U(r/n?7,)"@"-)F (r)sin[ (n6 —)/(n—1)] 


ow 1 o @ nf(2n—2) 
8.0) = 20(1— ayer F(R} 
, ce 0—n (2n—1)0—a 
u(r, 0) = P(rjoos(—7)—Q T (| 


u(r, 8) = Palrsin(7—2) + Q,()sin{ EE" 


where 
nD, (2) = (12) P(r /ry)O-2{(n— 1) 7 F(x) + 0Fy(0)} 
(2n—1)n™""-DQ,, (7) = —(1—r)-P(r/1)4O"-9(n— 1) F;,(7) ; 
and F,(r) = F{A, B; (2n—1)/(n—1); 7}, 
A+B =[n/(n—1)]—B, AB= —Bn(n—})/(n—1)? 
provided that n > 2, a >0>n/n. 
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When @ = z, we have argz = 0 and ¢ = 0. Again, when 0 = a, 
we have argz = « and ¢ = 0. Consequently, the ‘corresponding’ 
compressible flow solution given by (8) does, in fact, represent the 
flow of a compressible fluid between the two walls argz = 0 and 
argz = « regarded as rigid boundaries. 

In the case when the walls are at right angles, i.e. when n = 2, 
the solution admits the . form, 


4ih(7, 0) = U(r/7,)F,(7)sin 26 
—44(7, 0) = U(r/7,) ‘on Pf F'4(r) +F,(7)}cos 20 


u(r, 0) = Q,(r)cos 36—P,(7)cos 6 

y(7, 0) = Q,(r)sin 30+P,(r)sin é . 

where 4P,(r) = (1—r)-P(z/7,) {or Fo(7) + 2F,(7)} 
—12Q,(r) = (l—2)-A(z/1y)!7 Fal) 

and (choosing y 4) 6B = 24,A = 24,B=—3 

and P(r) = 1—(5/2)r+ (35/16)2— (21/32)r°. J 


The equations from which the cusps are to be determined in this 





case are, using (4), 

(i) palblibep al xe = 0 

(ii) 7?(13,2307°—72,76574+- 160,96573— 178,7527?+- (10) 

4.99,5847—29,972) == @, 

The latter of these two equations has no root between 0 and 1 
(excepting tr = 0), and is ignored in consequence. The first equation, 
however, yields the root + = 0-31146 (approx.). In order to find 
the value of @ corresponding to this value of 7 we notice that equa- 
tion (10) (i) is mF, (7)+2rF),(7) = 0 with m= 2. Using (3) with 

m = —% = —$}m we find that, since 7 >0 > 47, 0= 3x. The 
value of the critical streamline is found, using (9), to be 
—0-032215(U/7,). 
Again, the maximum value of 7 on this streamline is easily verified 
to have the same value of 7 as corresponds to the cusp, viz. 
= 0-31146. 
Corresponding to this value of 7, the Mach number attains the 


maximum value of 1-504 (approx.), without violating the principles 
of continuous flow. 
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When the fluid is assumed incompressible, denote by (q;,6;) the 
values of (q,@) at the point (r, x) in polar coordinates. Then, using 
(6) with n = 2, 

(i) r= q,/2U, Gi) x= 7-4. © (11) 
Introducing x’ and 6;, where x’ = x—}m and 6; = 37—8,, we have 
%,=x' where |6;|<4}, |x’| <<}. (12) 

When the fluid is assumed compressible we have, using (9), 

(i) r? = P?+Q2—2P, Q, cos 40, 

(ii) tany = (Q,sin 30+ P,sin 6)/(Q, cos 30—P, cos @). (13) 
Introducing 6’, where 6’ = ?7—@ (|0’| < 37), we write the latter as 
P,sin Ad = Q,sin(3A0+4%), where A@ = 6’—6%. (14) 

This quantity A@, which represents the change in the direction of 
the fluid motion due to compressibility, is positive in the region 
ba < argz < }n, negative in }7 > argz > 0, and zero along the 
lines argz = 0, }x, and $7. Using (14), it can be seen that the curves 


of Aé@ plotted against x’ (= 6;) for constant values of 7 (excepting 
7 = 0), have maximum values A@,, given by 


cosec AO, = P,/Q. (7 #0). (15) 


Again, if x}, is the value of x’ corresponding to Aé,,,, then 
8xin = 7—6A8,,. (16) 


These values, correct to the nearest minute, are tabulated for various 
values of 7 up to tr = } in Table 1. 

Using (13), we find that, when x = }z, r = P,+ Q, = q,/2U, and, 
when xy = 0 or 47, r= R—Q, = 4q,/2U. For these values of x, 
I have illustrated in Table 2 the numerical relationship between 
vt = (q/@) and (q;,/q) for various values of +. These relationships 
form the inner and outer bounds of the deviation in the velocities 
due to compressibility for any value of x. 

It should be noted that, using the familiar principle of reflection, 
the relations (8) also represent the solution for the symmetrical flow 
between the two rectilinear walls y= +axtana (0 <a < }n), 
where the flow at infinity is directed towards the origin along y = 0 
in the sense of x decreasing. 

In conclusion I wish to take the opportunity of thanking Mr. H. 
Williams for his untiring work in checking the numerical results. 
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TABLE 1 





0-02 0-04 0-06 0-08 0-10 0-12 0-16 
30’ a 1° 36’ | 2° 16’ | 2° 59’ | 3° 46’ 5° 36’ 
22° 8’ | 21° 44’ | 21° 18’ | 20° 48’ | 20° 16’ | 19° 41’ 18° 18’ 






































TABLE 2 





+ 0-02 | 0-04 | 0-06 | 0-08 | 0-10 | 0-12 | 0-14 





=) 0-1414 | 0-2000 | 0-2449 | 0-2828 | 0-3162 | 0-3464 | 0-3742 
(qj/@)y<o || 0-1367 | 0-1864 | 0-2186 | 0-2433 | 0-2603 | 0-2719 | 0-2789 
(q;/%)y=4n\| 0°1390 | 0-1932 | 0-2324 | 0-2634 | 0-2888 | 0-3101 | 0-3280 



































REFERENCES 
1. G. Temple and J. Yarwood, M.A.P., ARC. Rep. and Mem. 2077 (1942). 
2. A. Chaplygin, Ann. Sci., Moscow Imp. Univ. (Math-Phys. Sec.) 21, 1. 











ON PARALLEL FIELDS OF PARTIALLY 
NULL VECTOR SPACES 


By A. G. WALKER (Sheffield) 


[Received 9 August 1948] 

1. Introduction 

In the present paper I discuss parallel vector spaces, or planes, in a 
tiemannian space V,, and concern myself particularly with planes 
which are partially null as defined in § 2, for these appear to have 
received very little attention. It may have been assumed that 
properties of non-null planes hold also, by a kind of continuity, 
for partially null planes, and that the latter do not therefore merit 
special attention. That this is not the case is seen at once when we 
observe that, whereas V, is reducible if there exists a parallel field of 
non-null vectors (or more generally non-null planes), the space is not 
necessarily reducible if the vectors of a parallel field are null.+ 
Another significant difference is that, whereas a parallel field of 
non-null planes implies the existence of a symmetric second-order 
tensor, other than the fundamental tensor, whose covariant deriva- 
tive vanishes, the same cannot be said of a parallel field of partially 
null planes. This and similar theorems will be proved in the present 
paper. I also define and discuss recurrent tensors and the vector spaces 
generated by them. Special cases of recurrent tensors are tensors 
whose covariant derivatives vanish. 

In a later paper I hope to give a canonical form for the metric of 
a Riemannian space admitting a parallel field of partially null planes 
of any dimensionality. 

The literature on subjects relating to parallel vector spaces (mostly 
non-null) is extensive, but attention may be drawn to papers by 
L. P. Eisenhart (1923)t and H. Levy (1926)§ on tensors whose 
covariant derivatives vanish, by L. P. Eisenhart (1925|| and 193877) 
on parallel vector fields, by T. Y. Thomas (1939){{ on the decomposi- 


+ This fact can be seen to follow from the canonical form given by L. P. 
Risenhart, Annals of Math. 39 (1938), 316, for the metric of a space admitting 
a number of independent parallel vector fields. It is also illustrated by the 
example given in § 4 of the present paper. 


t Trans. American Math. Soc. 25, 297. § Annals of Math. 27, 91. 
|| Trans. American Math. Soc. 27, 563. tt Loc. cit. 
tt Monats. fiir Math. Phys. 47, 388. 


[Quart. Journ. of Math. (Oxford), Vol. 20, Sept. 1949, pp. 135-45] 
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tion of Riemannian spaces in the large, and by Makoto Abe (1944)7 
on the relation between parallel planes and Cartan’s holonomic 
group. Mention should also be made of a paper by Yung-Chow Wong 
(1945) on quasi-orthogonal ennuples. This is concerned with par- 
tially null sets of vectors, and it is possible that quasi-orthogonal 
ennuples could be used with advantage in the study of parallel planes 
in general. 

There is no generally accepted terminology and notation in the 
subject of the present paper, and I therefore give precise definitions 
where they appear desirable. Certain elementary properties of planes 
are stated without proofs. Suffixes i, j, k take values 1, 2,..., n 
throughout; suffixes a, b, c are used to distinguish between different 
vectors of a system and take values according to the dimensions 
of the system. The fundamental tensor of V, is g,;, and we need 
only assume that these components possess first derivatives with 
respect to the coordinates 2‘ of an allowable system, and that 
ll Giz || A 9. 


2. Planes at a point 

An r-plane at a point P of V, is an r-dimensional algebraic vector 
space whose elements are vectors at P. The vectors of a plane are 
all covariant or all contravariant; from a covariant (contravariant) 
plane we obtain the equivalent contravariant (covariant) plane by 
raising (lowering) suffixes in the usual way by means of the funda- 
mental tensor. Any set of r independent vectors of an r-plane is a 
basis; if A* (a = 1,...,r) is such a basis, the r-plane is the set of all 

a 
vectors aA’. The dimensionality r of a plane satisfies 0 <r <n. 
a 


The only 0-plane is the zero vector, and the only n-plane is the 
tangent plane, i.e. the set of all vectors at P. 

Let p, p’ be two planes at P. Their intersection, pn p’, is the set 
of all vectors common to p and p’, and their sum, p+-p’, is the set 
of all vectors obtained by adding a vector of p and a vector of p’, 
the zero vector being regarded as belonging to every plane. Both 
pnp’ and p+ ’ are planes, and if the dimensionalities of p, p’, and 
pnp’ are r, 8s, t respectively, then p-+-p’ has dimensionality r+s—t. 

Two planes are orthogonal if every vector of one is orthogonal to 


+ Proc. Imp. Acad. Toyko 20, 56, and 177. 
t Annals of Math. 46, 158. 
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every vector of the other. The plane conjugate} to an r-plane p is 
the (n—r)-plane of all vectors orthogonal to p. If p’ is conjugate to p, 
then p is conjugate to p’, and p, p’ are orthogonal. The 0-plane and 
the (tangent) n-plane are conjugate. 

A plane is null if it is self-orthogonal. Thus every vector of a null 
plane is null, i.e. satisfies g,;A‘’ = 0, and is orthogonal to every 
other vector of the plane. Conversely, if the basis of a plane consists 
of mutually orthogonal null vectors, the plane is null. The 0-plane 
is null. 

The null-part of a plane p is the intersection of p and its conjugate, 
and is a null plane. The null-part of a null plane is the plane itself. 
A plane, other than the 0-plane, whose null-part is the 0-plane, is 
said to be non-null. The tangent plane is non-null. A plane which 
is not non-null is said to be partially null. Thus a null plane is 
partially null, and a partially null plane other than the 0-plane has 
a null-part of non-zero dimensionality. 

Let the null-part p* of an r-plane p have dimensionality s (s <r). 
Then the plane p’ conjugate to p has dimensionality n—r, and the sum 
p+p’ has dimensionality n—s; this sum is the plane conjugate to 
p*. Since the (n—r)-plane p’ contains the s-plane p*, we have 
n—r >s. Hence, 


THEOREM 2.1. The null-part of an r-plane in V,, cannot have dimen- 
sionality exceeding n—r. In particular, the dimensionality of a null 
plane cannot exceed 4n. 

A basis of an r-plane is normal if it consists of r mutually orthogonal 
vectors each of which is either null or a unit vector. Every basis of 
a null plane is normal. If the null-part of an r-plane p is an s-plane 
p*, every normal basis of p consists of s null and r—s non-null 
vectors, the s null vectors forming a basis of p*. 

Let 7';' be a tensor of any order at P. Associated with this tensor 
and with a particular suffix, say i, are planes whose vectors 2' satisfy 
T;; A4* = 0; such a plane can have any dimensionality from 0 to n—r, 
where r is the rank of the matrix (7';'), in which i denotes the row 
and all other suffixes denote columns. The associated plane of maxi- 
mum dimensionality n—r is unique and consists of all the vectors * 
which satisfy 7;'A* = 0. 

+ We cannot now use the more familiar term ‘orthogonal complement’ 
because, if p’ is conjugate to a partially null plane p, then p and p’ are not 
complementary in the sense that p+’ is the whole tangent plane. 
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It will be found more convenient to consider not the planes asso- 
ciated with a given tensor but their conjugates. In particular, we 
define the plane generated by T;;' (for the particular suffix) as the 
plane conjugate to the associated plane of maximum dimensionality. 
Thus the plane generated by 7’; has dimensionality equal to the rank 
of the matrix (7';'). This rank is at most n, in which case the plane 
generated by the tensor is the tangent plane. 

The plane generated by a vector 7; is the 1-plane of vectors «7;. 
A tensor of order greater than 1 generates a plane for each suffix; 
these planes are all identical if the tensor is symmetric or skew sym- 
metric in every pair of suffixes, or if it has other symmetry properties, 
such as those possessed by the curvature tensor R,,;;,, which enable 
every suffix to be brought to the same position. 

Corresponding to an r-plane p is an rth-order skew symmetric 
tensor, the multivector of the plane. This is defined in relation to a 
basis At (a = 1,..., 7) by 

a 
parte = BarstrNh,, Nir = || ro! 
a dy ta@ u 
and the only effect of a change of basis is that the tensor may be 
multiplied by a scalar. This multivector is of rank r, and clearly 
generates the r-plane p. 
The measure of an rth-order multivector is defined as 
1 ae 
—D; i p Rove r| - 
lp Pie | 
When the corresponding plane is non-null and the multivector is 
constructed from a normal basis, its measure is 1. When the corre- 
sponding plane is partially null, the measure of the multivector is 
zero. 

More generally, if a partially null r-plane whose multivector is 
p®~t has for its null-part an s-plane whose multivector is p®~*, then 
pp; _;,can be contracted r—s times but no more without vanishing, 
and pimivterrntep, oo g = pip, 5, (1) 


where ¢ is a non-zero scalar. 
3. Parallel fields of planes 


A field of r-planes is a set of r-planes, one at each point of V,, 
possessing bases whose vector components 2 are differentiable func- 
a 
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tions of position. It follows that the first covariant derivatives 2‘, 
a 
of the basis vectors exist, and that the absolute derivatives 


DD. ~~ & 
dea 2" de 
exist along any curve of differentiable arc. 
A field of planes is parallel if, for any two connected points P, Q, 
a vector in the plane at P is displaced into a vector in the plane at @ 
by parallel transport along any curve from P to Q. For a fixed 
vector at P, the parallel vector at @ depends in general upon the 
curve from P to Q. If, however, the vector at Q is the same for all 
curves, and if this is true for every vector at P and for every pair 
of points P, Q, we say that the field of planes is strictly parallel. 
Strictly parallel fields call for little special attention, because from 
the definition it follows that a strictly parallel field of r-planes exists 
when and only when there are r independent parallel vector fields, 
i.e. vectors A’ (a = 1.,...,7) satisfying A*;, = 0. 
a a 


The field of zero vectors and the field of tangent planes are trivial 
examples of parallel fields. It is clear from our definitions that, if there 
are two parallel fields of planes, then others are formed by taking the 
intersection and the sum of the planes at each point. Also, since ortho- 
gonality between vectors is preserved by parallel transport, it follows 
that, if the planes of two parallel fields are orthogonal (conjugate) at 
one point, then they are orthogonal (conjugate) at every point. We 
can thus talk of the ‘intersection’ and ‘sum’ of two parallel fields, 
and of ‘orthogonal’ and ‘conjugate’ parallel fields. Further, since a 
normal basis at one pointisdisplaced by parallel transport into a normal 
basis at any other point, we see that the planes of a parallel field are 
either all null or all non-null or all partially null to the same extent, 
and that their null-parts constitute a parallel field. Hence 


THEOREM 3.1. Associated with a given parallel field of planes are 
(i) the field of all conjugate planes, and (ii) the field of ali null-parts. 
Both of these are parallel fields. 

The field of null-parts is, of course, the zero field if the given 
planes are non-null. 

If an r-plane p contains an s-plane p’ (s <r), and if p” is the 
(n—s)-plane conjugate to p’, then p n p” has dimensionality t > r—s. 
If p’ is non-null, pnp” and p’ are disjoint (pn p”N p’ = 0) and 
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(p nN p”)+p’ has dimensionality t+-s. But this plane is contained in », 
whence ¢+-s <r and therefore t= r—s. Thus (pn p’)+p' = p. 
If now p and »’ are planes of two parallel fields, p” and pn p” are 
parallel fields, and we have 


THEOREM 3.2. If a parallel field of r-planes contains a parallel field 
of non-null s-planes (s <r), there is a parallel field of (r—s)-planes 
orthogonal to the s-planes, and the r-plane at each point is the sum of 
the s-plane and the (r—s)-plane at the point. 


It should be noted that a decomposition of this kind cannot 
necessarily be carried out when the sub-planes are partially null. 
A decomposition can be effected at each point, but not by a process 
necessarily leading to a parallel field. In particular, when r = n, we 
see that a parallel field of non-null s-planes has a complementary 
parallel field of (n—s)-planes, the conjugate field, but that a comple- 
mentary field of this kind may not exist when the s-planes are 
partially null. This is at bottom the reason why V,, is reducible when 
it admits a parallel field of non-null planes but not necessarily when 
the planes of a parallel field are partially null. 


THEOREM 3.3. The field of planes with basis X' (a = 1,..., 7) is parallel 


if and only if 4, = AleX (2) 


for some functions A®,.. 

To prove this at any point P, let ¢’ be a unit vector at P, let C be 
a curve through P in the direction ¢’, and let y; be parallel along C 
and orthogonal to the r-plane at P. Then, if the r-planes are parallel, 
pu, A* = 0 at all points of C, and on differentiating along C we get 

a 


KG rN ;, tk = @ 
a 


at P. This holds for all directions ¢‘ and all vectors 4; orthogonal to 
the r-plane at P. Relations of the form (2) follow at once. 

The sufficiency of the above conditions is easily established by 
proving that any vector initially orthogonal to one of the r-planes of 
basis \* satisfying (2) remains orthogonal to these planes when 

a 


displaced by parallel transport along any curve. 
4, Parallel 1-planes 


A field of 1-planes must not be confused with a vector field. Its 
basis is a vector field, and from (2) with r= 1 we see that the 
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1-planes aA‘ are parallel if the vector field ‘ satisfies 

Mp = Nay, (3) 
for some vector a,. Replacing 4‘ by ¢A‘, where ¢ is scalar, has the 
effect of adding a gradient to a,. If a, is itself a gradient, we can 
write a, = —¢,/¢, and we then see that the vector field ut = ¢d* 
satisfies u*,, = 0, i.e. is parallel. Conversely, of course, the 1-planes 
whose bases form a parallel vector field are parallel. 

We observe that a, is always a gradient when 2‘ is non-null; for 
then Gig A'Nay, = Gag AW = 40, (Gij A’), 

Thus a parallel field of non-null 1-planes implies the existence of a 
parallel vector field. 

When 2° is null, however, a, need not be a gradient, and there 
does not necessarily exist a parallel vector field. An interesting 
consequence is that when A! at P is displaced by parallel transport 
round a closed curve and returns to P, the final vector is ¢A*, where ¢ 
depends upon the path and is not in general unity. 

To prove that a, is not always a gradient when 2° is null, let V, 
have coordinates x, y,z and metric 


ds? = u(x, y, z) dx? +-dy?+ 2dadz 
and consider the null vector field A‘ = 8§. It can easily be verified 


that 9 
Mii, ya be 
2 dz 
Thus the 1-planes ad‘ are parallel, but a, is a gradient only if % has 
the special form zf(x)-+-g(2, y). 


5. Recurrent tensors 
A tensor 7';° of any order will be called a recurrent tensor if it 
satisfies AS aes 
Tin = Ti (4) 
for some vector a,, the vector of recurrence.t From (3) we see that 
the basis of a parallel field of 1-planes is a recurrent vector. A special 
+ The term ‘recurrent’ was first used in this sense by H. 8. Ruse (Proc. 
London Math. Soc., in press) in his study of spaces of recurrent curvature, 
i.e. spaces whose curvature tensors satisfy Rpijgz, = Rnijz, x1. Here x; turns out 


to be always a gradient, a property not possessed in general by the vector a, 
in (4). 
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case of a recurrent tensor is one whose covariant derivative vanishes, 
the vector a, then being zero. 

If 7';' is recurrent, then, for any scalar ¢, S; = ¢7';' is also 
recurrent. In particular, ¢ can be chosen so that S;', = 0 when 
and only when a,, is a gradient, ¢ then being given by a, = —¢,/¢. 

It is clear that every product of recurrent tensors is a recurrent 
tensor. Since the fundamental tensor is recurrent, it follows that 
every contraction of a recurrent tensor by means of the fundamental 
tensor is recurrent. 

I now give a number of theorems showing how recurrent tensors 
are related to parallel fields of planes. 

THEOREM 5.1. The field of planes generated by a recurrent tensor is 
parallel. 


It is sufficient to prove that the field of conjugate planes is parallel. 
Let 7';' be a recurrent tensor, and let A‘ be any vector satisfying 
T; Xi = 0. 
Differentiating covariantly and using (4), we have 
Tr, = 9. 
Since every solution of (5) is a vector of the conjugate plane, it 
follows from (6) that A‘, for each k is in the conjugate plane. The 


basis vectors of this plane therefore satisfy relations of the form (2), 
i.e. the conjugate planes are parallel. Hence the planes generated 
by 7’; are parallel. 

Conversely, it follows from the next theorem that every parallel 
field of planes is generated by some recurrent tensor. 

THEOREM 5.2. The multivector corresponding to a parallel field of 
planes is recurrent. 

This is easily verified from the definition of a multivector and 
from (2). We find 

pit , = peta, Uy, = Ary. (7) 

The above two theorems provide an alternative test for a parallel 
field of planes. We have 

THEOREM 5.3. A necessary and sufficient condition for a field of 
planes to be parallel is that the corresponding multivector should be 
recurrent. 

This condition is necessary because of Theorem 5.2, and sufficient 
because of Theorem 5.1. 
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Before leaving the subject of multivectors it may be of interest to 
observe that the vector a,, in (7) is effectively due to the null-parts 
of the parallel planes. For, if p*~* and p*~“ are the (recurrent) 
multivectors corresponding to a parallel field of planes and to the 
field of null-parts, and if a,, ag are their vectors of recurrence, then 
each side of the identity (1) is a recurrent tensor, the vectors of 
recurrence being 2a, on the left and 2af+-¢,/¢ on the right. Hence 
d),—a@j, is a gradient, and we have 

THEOREM 5.4. The multivectors corresponding to a parallel field of 
planes and to its sub-field of null-parts can be chosen so that they have 
the same vectors of recurrence. In particular, the multivector corre- 
sponding to a parallel field of non-null planes can be chosen so that its 
covariant derivative vanishes. 


The latter case is when the multivector has constant measure. 


6. Recurrent second-order symmetric tensors 

We finally consider tensors 7; which satisfy 7, = 7, and 
T;; ;, = T,;a,. There is always a tensor of this kind with a, = 0 
associated with a parallel field of non-null planes, and the following 
theorem is well known. 


- THEOREM 6.1. Let X* be a normal basis of a parallel field of non-null 
a r ‘ 
r-planes, and write T;; = Ty, = > egdA,A;, where eg = gyrA = +1. 
a=1 aa aa 
Then T,; generates the given r-planes and satisfies T;; ,, = 0. 

This can easily be deduced from (2). It also follows from the fact 
that, if pr is the multivector of unit measure corresponding to the 
planes, its covariant derivative is zero, and 7;,; is given by 

l i a5 
Ty = Ga Pe Patent 

Turning now to parallel fields of partially null planes, it is natural 
to seek for such a field a theorem analogous to Theorem 6.1. It 
appears, however, that no such theorem exists, even when the 
requirement is weakened to that of a recurrent tensor. I prove 


THEOREM 6.2. There is a parallel field of partially null planes which 
is not generated by any symmetric recurrent tensor of the second order. 


Consider the V, with coordinates x, y,z and metric 
ds? = y? dz*+dy?+ 2dzdz, 








144 A. G. WALKER 


a special case of the V, in §4. The null vectors ad{ are known to 
form a parallel field of 1-planes, and there is therefore a parallel field 
of partially null 2-planes, the conjugate field. A tensor 7; = 7;; is 
associated with these 2-planes if and only if 7;, = 0 (¢ = 1, 2, 3) 
and rank (7;;) = 2. We therefore seek a tensor 7;; satisfying 
T= Tj, Te =9, TMaTe—(TrP? #0, Typ = Tyga, 
(i,j, = 1, 2,3) (8) 
for some 4,,. 

Suppose firstly that a, is a gradient, say —¢,/¢. Then by consider- 
ing $7;,; in place of T;; we have an equivalent problem with a, = 0. 
We can therefore assume a, = 0 in (8). The only non-zero Christoffel 
symbols are found to be 


ae (3 
pease MS ha 


and the equations 7;,;;, = 0 become 6; Th. = 0 (k = 1, 2, 3), and 


0, Ty = —2yT 2, 0, Ty. = —yT x2, 
02 Ty = 03T yy = 02 Tyg = O5Ti2 = 9. 
It is at once seen that these equations have no non-zero solution. 

Suppose now that a, is not a gradient. We find 0,7), = Th. a,, 
whence 7}, = 0 since a, is not a gradient. We now have 6, 7), = 7}. a,, 
and again 7,, = 0, which contradicts 7, Th. ~ (Ty)? There is there- 
fore no solution, and the theorem is proved. 

We thus see that parallel fields of planes fall into three classes: 
(i) those generated by a second-order symmetric tensor whose 
covariant derivative vanishes, (ii) those generated by a second-order 
symmetric recurrent tensor whose vector of recurrence is not a 
gradient, and (iii) those which are not generated by any second- 
order symmetric recurrent tensor. All three kinds exist, but except 
for parallel fields of non-null planes, which belong to class (i), it is not 
clear what properties, other than those used to define the classes, 
distinguish between parallel fields of different classes. 


7. Recurrent second-order symmetric tensors (continued) 

We know from Theorem 5.1 that, if 7;; is recurrent, it generates 
a parallel field of planes. This field is, however, non-trivial (i.e. is 
not the field of 0-vectors) when and only when 7;; is singular, 
ie. || Z;; || = 0. I shall now prove that there exists a parallel field of 
planes even when 7;; is non-singular. 
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THEOREM 7.1. If there is a recurrent tensor T,; other than a multiple 
Of 943, then there is a non-trivial parallel field of planes. 


The case when ||7;,;|| = 0 is covered by Theorem 5.1. Suppose 
now that || 7;;|| A 0 and that 7); = T;;a,, and write ¢ = || 7; |I/\| 94; |l- 
Then ¢ is a non-zero scalar product of T’s and g’s of degree n in the 
T’s, whence ¢;, = nda,, and a, is a gradient. Writing S,, = ¢-¥"7;,, 
then S;;, = 0. We have thus proved the existence of a tensor S,; 
whose covariant derivative vanishes; S;; is non-singular and is not 
proportional to g;;. 

Consider now the polynomial 


Each coefficient 7, is a sum of scalar products of S’s and g’s and is 
therefore constant because S;;, = 0. Every root of the equation 
7(p) = 0 is therefore a constant, and is non-zero because |{§;;|| A 0. 
If p is any one such root, write 
U5 = p93 —Sy- 

Then U;;;, = 0 and || U;; || = 0. The field of planes associated with U;; 
is therefore parallel and non-trivial, and our theorem is established. 

The above theorem suggests that in a search for tensors Tj; satis- 


fying 7;;,, = 0, or more generally 7; ;, = J;,;a,, we should first seek 
all the parallel fields of planes; if there is no such field other than 
the two trivial fields, then there is no recurrent second-order sym- 
metric tensor other than g;;. 


3695.20 
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1. Introduction 

Ir A is any matrix, we denote by A’ the transposed matrix and by 
A the conjugate complex matrix, while |.A| denotes the determinant 
and o(A) the trace of a square matrix A. By [q,...,a@,] we designate 
the diagonal matrix with diagonal elements 4,,...,a, -and all other 
elements zero. Also J denotes a unit matrix and O a zero matrix 
whose order will be clear from the context. 

If X = (x;,) is a real nxn matrix, then the n* coordinates 
24, (j,k = 1,...,n) determine a point in the Euclidean space Sp of n? 
dimensions. The manifold of Sp defined by X’X = I is of dimension 
in(n—1). We call this manifold VM; it is the space of orthogonal 
matrices. Let v(D) be the Euclidean volume of the domain D; then 
we shall give an alternative proof of Hurwitz’s (1) evaluation* of 
v(Mo). 


n ir 

Ls 7 
THEOREM 1. (Mo) = 2448) | | Tar)’ 
r=1 . 


Similarly, if Z = (xj,+-ty;,) is a complex nxn matrix, the 2n? 
coordinates 2%, y;, (j,4 = 1,...,%) determine a point in the Euclidean 
space S, of 2n? dimensions. The manifold of S, defined by Z'Z = I 
is of dimension n?. We call this manifold M; it is the space of all 
unitary matrices of order n. Let the sub-manifold of M, defined by 
|\Z| = 1 be called MY. It is of dimension n?—1. Hurwitz (1) 
evaluated v(M/{?), and his method may be modified to give v(My). 

By our method we find 

n 

THEOREM 2. v(My) = 2inin+0 


r=] 


TT 


Tr) 
Siegel (7, 8, 9), using a different definition of volume, has also 
discussed the volume of orthogonal and unitary space. 


* Actually Hurwitz finds the volume of the space of proper orthogonal 
matrices. 


(Quart. Journ. of Math. (Oxford), Vol. 20, Sept. 1949, pp. 146-54] 
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2. Orthogonal space 

Let (r,,72) be the rth pair of integers in the sequence (1, 2), (1, 3),..., 
(1, ”), (2, 3),..., (2,2),..., (n—2,n—1), (n—2,n), (n—I1,m), and (8,, 8) 
the sth pair. Let the elements of the nth order skew-symmetric 
matrix X above the diagonal be re-indexed as the vector x of 4n(n—1) 
We define the symmetric matrix C of order 


rows with 2, = 2,,,,- 


yn(n—1) by 2’ Ox = o(A'X'AX). 


Lemma 1. If A is symmetric, then 
6 — 2A®), |C| ee 2inin—-1)| A |”-1, 
where A® (or C,(A)) is the second compound or adjugate of A. 


Proof. oA’X'AX) = 333) 4.2... Go lnm 


i <2 
cul. - > p vp, ra(Fr, 81 Dy. 59 Ay, 5,4 r5s,)%o.85 
T1<T2 81<82 


ya Le v, Chg Ug, 


where c.. = 2 Fy 9 D8, =C¢ 
rs a ‘sr? 
Aro, r282 


since A’= A. That is: C = 2A®., Also A@ is of order $n(n—1) 
and* | A®@| = |A|"-1, 


i {aT} 
LEMMA 3. Let E,,(«) = | orb 
Ln 
where T is real skew-symmetric of order n, ,, is the real Euclidean space 
of 4n(n—1) dimensions with coordinates 
tp (l<j<k<n) and {dT} = TT dtye. 
i< 
Then, for « > —1}, 


ie Ee weet 


I+7, t 
—~ 1)’ 





Proof. Putting I+T = ( 


we have 


\I+T| = I+T|\1+¢U+7)%), {dT} = Tay}, 
where {dt} = TL dt;,. Then, if v,_, is the space of (n—1) dimensions 
j=1 


* See Wedderburn (3) 66 or Macduffee (5). 
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with coordinates t;,, (j = 1,...,n—1), we have 


Fe {at} {ay} 
gc | ( J ORs) ee 


Mn-1 Vn-1 





Now put ¢ = (J+ 7,)u; then 
U(I+-7,)-% = w'(I—T,)u = u'u 
since 7 is skew-symmetric. Also 
{dt} = |I-+T,| {du} 
{du} ah(n-1) T[a+ b(m-+- +$(n-+1)] 1 
and loathe - lat+n) ” (1) 


Vn-1 


provided that « > —}(n+1). Hence 


- a aoe et) 1)] {dT, 
E,,(«) = ait ee - (J+ = 4 waded -l1 


a —t- n) 


= on Tlet Het], E,,-4(«). 








T(a+n) 
Also, if « > —14, 
(__ae P(a+$ 
E.(«) = : — 2) 
- re * T(a+2) 


Hence the result. In . 


atk ne —+t) 9—Hn—2Xn—1) 
B(—) = [7 Pop a aa 1 I rear 


Proof of Theorem 1. We arrange the 2” (= N) diagonal matrices 
[+1,+1.,...,-1] in a sequence 
J, = [1,1, 1... 1], 
J, = [—1,1,1,..., 1], 
J, = -~ wal nae 
J, = [1, —1, 1...., 1] 


and so on, in such a way that a is for med from J,_, by changing the 
sign of only one diagonal element. Since* 


N 
> ,X+1| =N, 
r=1 


* This result is due to Lipschitz; see Muir (10) 289. 
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for each X there exists at least one J, such that 
\J,X+1| #0. 
We define sub-manifolds of Mp thus: M, is the set of points X of Mp 
such that 
|\J,X+1I|=0 for p=1,2,...,.r—1, but |J,.X+J| #0. (2) 
Then M,, M,,..., My are non-overlapping sets and 
My = M,+M,+...4+My, 
so that v(Mo) = v(M,)+v(M,)+...+0( My). 
If X belongs to M,, put 
T = (I+J,X)*(I—J,X); 
then, as X’X = I, we have the Cayley transformation 
X = J(I—T)I+T), 
where 7” = —T and, when r > 1, by (2) and (4), 
\J,F(L—T)(I+T)*+1| =90 (p = 1,2,....7—]), 
i.e. \(I+-J,4,)+(1—J,d,)T| = 9 (p = 1,2,...,.7—1). (5) 
First we evaluate v(M/,). For all X in MU, 
X(I+T) = I-T, 
so that, if dX, dT’ denote matrices of the differentials of the elements 
of X, T respectively, we have 
dX (I+-T)+X dT = —dT, 
dX (I+ T) = —(1+-X) dT = —2(1+T)71 dT. 
Hence dX = —2(14+T)71 dT (I+ T)-1. 
Now the Euclidean metric is given by 
ds* = o(dX'dX) 
= 4o{(I+ 7") dT’ (I+7")7(I+T)+1 dT I+T)7} 
= 4o{([+7T’)+ dT’ (I+TT")“ dT} 
= dt’ G dt, (7) 


where G is symmetric and dt denotes the vector formed by rearranging 
the elements of the skew-symmetric matrix dT’ as in Lemma 1. Hence 
the Euclidean volume element [(6), 60] for the manifold M, is 


dv = |G|* {dt}, (8) 
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where, by Lemma I, 
|G| = ginin—D| J+ TT’ |i a peed s 
+75 
{aT} 
[+7 h> 
Hn 
= Qin) FP (—1) = 2Hn—1kn+4) 1 re 


v(M,) = Qin(n-l) 


ait 
T(4r)’ 
by Lemma 2 
Since \(I+J,)-+(1—J,)7| = 219° 4|=0, 
Ol 
we have the same domain of integration p,, and volume element dv 
for v(M,), so that v(M,) = v(M,). (10) 
If r > 2, some of the equations (5) are identities, but, since 
I gdh wn TM Ny 8, Bing ia Bd 
with exactly two negative elements in the jth and kth columns (say) 
\(1+J,-24,) + (J.J) P| = 2th. 
Thus M, is mapped by (4) into a manifold m, of the space defined by 
(5), and these equations include the hyperplane ¢;, = 0, so that m, 
is at most of dimension 4n(n—1)—1. Also we get the same volume 
element for MV, as in (8). Hence, when r > 2, 
»(M,) = | IG|t {dT} = 0. (11) 


ty 


The result follows from (3), (9), (10), (11). 


3. Unitary space 
Lemma 3. If A, B are real matrices and 
C=A+iB, |A| £0, 

2 


then* BA 


|= 1601. 


Proof. 
A -B\ (I O/A O I 
BA ) ~ \BA I}\O A+BA-—1B]\O 


* Cf. Murnaghan (11) 62. 








Hence 
|A —B 
IB A 





|CC\, 


since 
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|A||A+BAB| = |A||(A++B)(1—tA“*B)| 


\C| = |A||I—éA-1B]. 


Lemma 4. Let u, v, x, y be real n-rowed vectors, the matrices A, B, C 


as in Lemma 3, and 


w= ut+iv, z= 2+, a’ = (u'v’), v = (s’9}), 
then, if w= Cz, 
we have (81, S29---» San) _ OG), 


O46, s Sigcenp as 


Proof. We have 


u = Ax— By, v= Ba+Ay, 
i.e s= A -B t 
.e. “ts or 
ee a A -Bl was 
Thus gh at et Dt = |CC\|. 
" x .43° et 








Lema 5. If f(A) is a polynomial in the elements of the matrices A 
and A such that f(AB) = f(A)f(B), then either f(A) is identically zero, 
or f(A) = AA", where A = |A|, and g, h are non-negative integers. 

This is a known result [Murnaghan (11) 33, 194, or Weyl (12) 26, 
264). 

Lemma 6. If K, is a skew-Hermitian matrix, we denote by k, the 
real vector whose n® components are the real and imaginary parts of the 
elements of K, above and on the principal diagonal in some fixed order. 
If A is any n-th order matrix, we define a matrix C = C(A) of order n? 
by o(A’R’, AK,) = ki, Che. 

Then C is a real matrix of determinant 2%"-)|AA|", and is symmetric 
if A is Hermitian.* 

Proof. Firstly, by direct calculation, we have 

o(K; K,) = ki Jk, 
where J is a fixed diagonal matrix containing in its diagonal n 1’s 
and n(n—1) 2’s, so that |J| = 2%"-), Next, if K, is skew-Hermitian, 


* We thank the referee for his modification of Lemmas 5 and 6. 
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then AK, A’ is also skew-Hermitian, and, since the correspondence 
K, - AK, A’ is obviously linear, if we put AK, A’ = Ky, we have 

ks = F(A)kg, 
where F(A) is a real n?-rowed matrix whose elements are bilinear 
forms in the elements of A and A. Also F(AB) = F(A)F(B), so that 
| F(AB)| = |F(A)||F(B)|. Now F(J) = J, so that |F(A)| is not 
identically zero. Hence, by Lemma 5, |/'(A)| is a product of integral 


powers of |A| and |A|, and a senignelal of degrees shows that 
|F(A)| = |AA|". But 


kCk, = o(A’ K, AK,) = of K,(AK, A’)} = k Jkg = ky, JF(A)heo. 
Hence C= JF(A), a real matrix of determinant 2”-)|4 A|*. 
Finally, if A is Hermitian, A’ = A; so that 
ki, Ck, = o( A’ K;, AK,) = o( Ky A’K, A) 
= o(A’K, AK,) = k Ck, = K,C’k,, 


" 


and se CU = ¢ 
{dA\fdS' 
Lemma 7. Let F(a) = | —+ > 
‘ |I+ KK’ |=’ 
on 
where K is skew-Hermitian, K = A-+-iS (A, S real of order n), o,, is 
the real Euclidean space of n* dimensions with coordinates 





a, (l<j<k<n), &, (l<jck<cn) 
and {dA} = [J da,,, {dS} = [J ds,,.. 
j<k I<k 
Then, for « > —}, 
F (x) — 2-n(n+2a—1) aT (2 a+r) 3 
4 T (a--r)}? | 





Proof. Putting 


I+K (Fe ae } K, = A,+i8,, k =a+is, 
we have \I+K| = |J+ K,|(1+«+is,) 
where Kk = k'(I+ K,)“1k. 
Also {dA} = {dA,}{da}, {dS} = {d8,}{ds} ds, 


n—1 fo 
ran ‘da} = Ci al , 
where {da} = | l da;,; {ds} = | [ ds ;,. 
, en 
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Then, if r,,_, is the space of 2n—2 dimensions with coordinates aj, 8jn 
(j _ I josey n—1l), we have 


pes ad (44,}{aS;}, 


= J oe 


— 


where 








_ wD (a+n—}4) {da}{ds} 
~ Ta+n) (1 x)2+2n—1? 


Tn-1 


ifa > 4—n, Put 

k = (I+ K,)w, w= ut+w; 
then k= 0 ([—K,)w = Ww = w'u+v'v 
since K, is skew-Hermitian. Also, by Lemma 4, 


{da}{ds} = |I+ K, Kj| {du}{dv} 
so that, by (1), 








(K,) = mT(a+n—})|I+K, Kj! | {du\{dv} 
baie ie se P'(a+n) (1+-w’u+v'v)+2n-1 
a an Tot e—g ate) 7 14K, Ki) 


T(a+n)P(2a+2n—1) 
provided that « > —4n. Thus 
pale Ns a | An 








F, (a) = Tatn)r (2a+2n—1) n—1(%) 
221 —a—M apn D'(2a-+n)F,, ~1() 
* Sees. {['(a+n)}* 
Also, if « > —4, 0 
za ds = ml(2 x+ 1) 
Ala) = | Oet — Path 


= 


The result now follows. In particular 


F,(0) — 9-n(n-1) I] mm 
r=1 


Proof of Theorem 2. The proof is similar to that of Theorem 1. 
We consider sub-manifolds L, of My, where L, is the set of points Z 
of My such that 


\Z+wl|=0 for p=0,l,...,(r—1), but |Z+w'l| 40, (12) 
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where w = exp{2zi/(n+1)}. Bochner (2) shows that for every Z of 
M, there is an r such that |Z+w"/| #0. We have 

v(My) = v( Lp) +0(L,)+0(L,)+...+0(L,). (13) 
If Z belongs to L,, we have the Cayley transformation 

Z = w'(I—K)(I+K) (14) 

where K’ = —K, and, when r > 0, by (12) and (14) we have 

\iI cotka/(n+1)+K|=0 (k= 1,2....,7). (15) 
As in (6) and (7), the Euclidean metric is given by ds? = dk’G dk, 
where dk is formed from dK as in Lemma 6, and also, by Lemma 6, 

|G] = Qn@n-)| 74 KK’ |-2, 


Hence, by Lemma 7, 


v( Lp) = | |G a A}ds} = 2mm | | (Fa): (16) 


ri 
On 


If r > 0, L, is mapped into a manifold J, of o,, defined by (15), so 
that J, is at most of n?—1 dimensions. Hence, if r > 0, 

o(L,) = 0. 
The result follows from (13), (16), (17). 
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PERTURBED FUNCTIONAL EQUATIONS 
By E. M. WRIGHT (Aberdeen) 
[Received 20 November 1948] 


1. Introduction 
In the theory of differential equations there are several results which 
connect the behaviour of the solutions of the linear differential 


equation* ¥ a, y(x) = 0, (1.1) 
in which x is real and the a, are independent of x, with that of the 
solutions of a perturbed equation 

Zz a, y”(x) — d(x, y, VT sors y"-»). (1.2) 
For example, if every solution y, of (1.1) tends to zero as x +00, 
if ¢ is a continuous function of x and of the y” which is 


7) (max ly!) 


ven 
uniformly in 2 for small y”, and if the initial values y”(0) are suffi- 
ciently small, then y(x) > 0 as x > +00. Thus any solution of (1.2) 
which is initially small enough remains small or, to use slightly 
different phraseology, the zero solution y = 0 of (1.2) is stable. In 
fact, more is true. If every solution y, of (1.1) tends to zero as 
« —> +00, it follows that every root of the algebraic equation 
F(s) = da,s” = 0 (1.3) 

has its real part negative, and conversely. If every such real part is 
less than —c, for some positive c, and if the other conditions are 
satisfied, we have 
ly(x)| = O(e*) (v <n) 
as «> -+oo. Poincaré (13), Liapounoff (10), Bohl (4), Cotton (7), 
Perron (11), and Bellman (1) have proved results of this type, some- 
times under rather stricter conditions on ¢. 

It is possible to prove a similar result for a wide class of functional 
equations. To illustrate the method I here consider the problem 
in which (1.1) is replaced by 


a » 
A{y(x)} =X [ ye—é) dk,(é) = 0, (1.4) 
wien 


* Here and throughout 7 is a positive whole number. Except where a 
further restriction is explicitly stated, v takes all integral values for which 
Oxsrven. 


[Quart. Journ. of Math. (Oxford), Vol. 20, Sept. 1949, pp. 155-65] 
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where n > 0,6 > 0, 


zx b 
k,(x) = h,(a)+ f i) dé, [ |5,{€)| dé <0, 
0 0 


each h,(x) is a step-function with at most a finite number of dis- 
continuities, and h,(x) has a finite discontinuity at x = 0. The 
function ¢ is replaced by a functional Q(x, y) depending on x and 
the values of ya) (0<E <b, v<n). 
My result concerns the equation 

A(y) = Q(@, y). (1.5) 
I select (1.4) and (1.5) to illustrate the method because they include 
as special cases the difference-differential equations in which I am 
primarily interested (see, for example, 16) and also because I am 
able to shorten my work considerably by quoting certain results due 
to Pitt (12) about the solutions of the equation 

A(y) = v(z), (1.6) 

where v(2) is a known function. 

The stability result for differential equations is usually proved in 
three stages. The first trivial stage replaces each nth-order differentia 
equation by a set of n simultaneous first-order differential equations. 
This device is inapplicable to my equation, except for a very special 
case, viz. the difference-differential equations in which 


A(y) = y™(x)+ 2. Pa y(x—b,,), (1.7) 


so that only one nth-order deriv wo e occurs.* It is, of course, 
possible to handle sets of simultaneous equations of the type of (1.4) 
or of (1.5) but, in’ the general case, there is no obvious gain in 
simplicity by replacing (1.4) and (1.5) by any sets of apparently 
simpler equations. 

The second stage in the differential-equation theory replaces the 
differential equation by an integral equation (or rather the set of 

* While this article was in preparation Dr. R. Bellman very kindly sent me 
a copy of a paper (2) awaiting publication in the Annals of Mathematics. In 
this he points out that (1.7) for general n can be replaced by n equations, each 
of the type of (1.7) with n = 1. He proves the special case of my Theorem 1 
for the equation 


y' (x+1) = ay(x)+by(x+1)- a puilylo a) }*{y(x+- 1)}', 


and the corresponding case of my gee om :. 
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first-order differential equations by a matrix integral equation). For 
my equation this is a much more substantial step, but it can be 
deduced very simply from Pitt’s elegant results (12). 

At the third stage, the existence, uniqueness, and stability results 
are deduced from the properties of the integral equation. In the 
differential-equation theory, the existence and uniqueness results 
are classical, but I have to consider them briefly for my equation 
(1.5). I prove the stability result by a very simple continuity argu- 
ment, which is also applicable to the simple case of the differential 
equation. 

It may be remarked that the polynomial F(s) of (1.3) is now 


replaced by the integral function 
b 


K(s) = > 8’K,(s) = > 8” J e-€ dk, (E). 

van vcn 0 
The zeros of K(s) are infinite in number, and it is naturally less easy 
to determine in any particular case whether all their real parts are 
less than some fixed negative number. For certain special examples 
of the difference-differential equation, this problem has been solved 
by Hartree and others (5, 8) and for other examples by Hayes and 
the author. 

I have also used the methods of the present paper to extend to 
equations (1.4) and (1.5) results proved by Cesari (6), Levinson (9), 
Weyl] (14), Wintner (15), and Bellman (1) about the boundedness 
of solutions of (1.2) when those of (1.1) are all bounded and ¢ 
satisfies appropriate conditions. In this case, the zeros of F(s) must 
all lie on or to the left of the imaginary axis, while the multiple zeros 
all lie to the left of this line. The appropriate generalization of this 
condition to K(s) is complicated by the fact that K(s) may have an 
infinity of zeros on or immediately to the left of R(s) = 0. For the 
equation (1.7), however, K(s) has only a finite number of zeros in any 
strip of finite width parallel to the imaginary axis, and the extension 
is a little simpler. I intend to discuss these results elsewhere. 


2. Notation 
In both (1.4) and (1.5) we may conveniently take ¥(x) to be the 
unknown function and suppose yx) defined by 


ya) = y(0)+ [ YE) dE (v <n). (2.1) 
0 
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We suppose the ‘initial conditions’ to determine the values of 
ye) (—b<a<0), y0) <n), (2.2) 
subject to the proviso that y™(x) is integrable Lebesgue in the 
interval —b < x < 0. If these values are such as to make 
ly(x)|< 8 (-b<xx<d0,v<n; —b<gx<0, v=n), 

we say that y satisfies a 5-set of initial conditions. 

In what follows s is a complex number and og is the real part of s. 
C always denotes a positive number, not always the same at each 
occurrence, independent of 

x, xX, 8, é, Y; 5, €, U, 

but possibly depending on any other parameters present. C,, C,,... 
are numbers of the type C, but C,, for example, retains the same 
value throughout. The constant implied in the O( ) notation is of 
the type of C. Both « and 6 are positive numbers to be thought 
of as small; « is supposed arbitrarily assigned, while the choice of 5 
depends on the value of e. 


3. Pitt’s results 

We begin by stating certain results about K(s) and the solutions 
of (1.6) which are (apart from trivial extensions and re-wording) due 
to Pitt (12). He studies a somewhat more general equation. Here 
we suppose k,(€) constant when € < 0 and when é > 6 and remember 
that k,,(€) has a finite discontinuity at é = 0. 


Lemma 1. There is ao, such that K(s) has no zeros for which o > o,. 


If o a Oo _ O1; then |K(s)| > C+ C\s"|. 


Pitt only proves the second part for o, < o < C, but the extension 
is trivial in our case. 
Lemma 2. There is a function w(x) such that, when o > og > 04, 
w 
{K(s)} = | es dw(a). 
0 
The functions w(x) exist and are absolutely continuous and w(0) = 0 
for0 <v < n—1, w(x) is of bounded variation in every finite interval 
and all the integrals 
[ e-e|dw(x) (v < n) 


0 
converge. 
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We now introduce a function (x) with the following properties. 
When —b < x < 0, f(x) = y(x) as defined by the initial values (2.2). 
When x > b, &(x) = Oforv <n. For 0 < x < b, the function (2) 
is chosen so that (x) is absolutely continuous for vy << n,0 <<a <b 


and so that \W(ar) | < Cs (v < nN; rare * “4 x < b), (3.1) 


whenever y satisfies a 5-set of initial conditions. We write f(x) = 0 
(~ < 0) and 


b 
fiz) = — > | (e—u) dk,(u) (x > 0), 
0 


so that f(x) = 0 for x > 26 and, whenever (3.1) is true, 
lf(z)| << C3 («x >0). (3.2) 
By Theorem 9 of the second part of (12), the solution of (1.6) with 
the given initial conditions is 


zx 


y(x) = Y(x)+ [ {fle—£)+-v(a—€)} dw(é) 
0 


x 


= Wla)+ | {f(€)+v(6)} deo(e—£) 


0 
for « > 0. Now w™(0) = w(0+)=0 for »<n. In general 
wr(0-+) £0, but we may put w”(0) = 0 without violating any of 
the properties of w(x) stated in Lemma 2. On this basis we see that 


x 


Ya) = Yax)+ | {FE)+-v(€)} dw(a—£) (vw <n), 
from which we have ' 


Lemma 3. The solution yo of (1.4) with the given initial conditions 
is given by 
F x 


yf (x “J fe) dwa—€) = | flw—€) dew™E) 


0 
and the solution y of (1.6) with the same initial conditions by 


x 


yx) —y p(x =f ve ) dw”(a—€) = [ w(a—€) dw”) 
0 0 
forally <nand all x > 0. 
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We may remark that y = w(x) is the solution of (1.6) for which 
v(x) = 0 (x < 0), v(z) = 1 (@ > O) and yY(xz) = 0 WY <n, x < 0). 


Also, for og > 0 andy < n, 
tio 1 d 
] 8”—*e as 
LSy,(v)(» AV) (> = 
Lw(a2+ )+w(a—)} = — ——— 
ai 1 Qari K(s) ’ 
Go —in 


but we do not use these facts here. 


4. Existence and uniqueness of solution 
We prove first an existence result for the solution of a fairly general 


functional equation. 


Lemma 4. Let 0 <a < Q,, let w(x) be an integrable function of x 
and let \w(a)| < C,. Let R denote the set of all continuous functions u(x) 


for which u(O) = 0, \w(ax)| < C3. 


Let x(x, u) be a function of x and a functional of u such that, for all 
U, Uy, Uy belonging to R, x(x, u) ts integrable with respect to x, 


Ix(w,u)| < Cy, 
and x(@, Uy)— xX, Ue)| < C; — \w(€)—w.(€)|. 
0<f<2r 


Then the equation 


u(x) = | x(€, w)w(a—€) dé 
0 


has a unique solution u(x) for 0 < x < GQ. 

The condition (4.1) is not required for the existence of a solution 
but only for its uniqueness. Let us write { for the right-hand side 
of (4.1). Then, provided that w(x) is of bounded variation, we may 
replace ¢ in (4.1) by any p(¢) such that p(f) > 0 as + 0. The power- 
ful ‘fixed-point’ methods of Birkhoff and Kellogg (3) will then enable 
us to prove the existence of u(x), but not its uniqueness. The 
corresponding relaxation of condition (5.1) may be made in Theorem 1 
and Lemma 5. In most cases likely to arise, however, (4.1) will be 
satisfied and we therefore assume this, since it enables us to use the 
simple device of successive approximations. 

We write u,(x) = 0 and 


uy(a) = | x(E, wys)ol(w—€) dé. 
0 
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We have at once \w,(a)| < O,C,2, 


|uy(x)—wya(2)| < Oy Oy f max juy_y(q)—ty2(n)] a, 
; 0<n<é 


N 1 2 
< CFC, Cy- Wr’ 
and so uy(x) > a limit u(x). 

If u’, wu” are two different solutions, let 


M(x) = max |u'(£)—u"(£)|. 
0<é<z 
We have - 
|u’(x)—u"(w)| < O,€, [ M(E) dé < C,0,2M(a). 
i) 
Hence M(x) < C,C;2M(zx), and this is impossible for 
2<i, < Cz'C5}, 
unless M(x) = 0. 


5. The stability result 
Our main result is 
THEOREM 1. If (i) o, < —ce < 0in Lemma 2 and if (ii) for some b, 
satisfying 0 <b, < b, for all x > 0 and for all functions y, y,, y, with 
absolutely continuous (n—1)th derivatives satisfying 
Ya)| <Cy YP) < Cy, ye) <O, (@ > —b, ¥ <n), 
Q(x, y) is a continuous, integrable function of x satisfying 
Q(x, y) = o(a) ( = aff) = max ye—8)) (5.1) 
se 
uniformly in x as «—> 0 and 
Q(x, y)—Q(@, y2)| << Cg max  |yf(E)—y3(E)|, (5.2) 
O<2=f<bi 
then, for every positive e, there exists some 5 = 5(e) > 0 such that every 
5-set of initial conditions determines a unique solution y of (1.5) which 
has an absolutely continuous (n—1)th derivative and also satisfies 
lyY(x)| << ee (v< n) (5.3) 
foralla > —b. 
We first prove 
Luma 5. If 0 < « < C, if the conditions of Theorem | are satisfied, 
and if (5.3) is true for —b <x <0 and all v < n, then yx) exists 


and is uniquely defined for0 <% <Q, v <n. 
3695.20 M 
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If we write y,(x) for the solution of (1.4) which coincides with y(x) 
for —b < x < 0 (except perhaps for x = 0, vy = n), we have 


Ya) = yP(a)+ | QE,y) dw) (w@>0; <n) (5.4) 
0 


by Lemma 3. If we take v = n—1 in this and write 


u(x) = y-Y(x)—yf"-M(e) and w(x) = w(x), 

z 
we have u(x) = | x(€, u)w(x—&) dé, 

0 
where w and x satisfy the conditions of Lemma 4 with C, = ),. 
Hence u(x) exists and is uniquely defined for 0 < « < Q. The same 
is true for y{"-) by Lemma 3 and so for y”“-». Hence, by (2.1) for 
v <n and by (5.4) for v = n, y(x) exists and is uniquely defined 
for x < G. 

Now, by Lemma 3, 


yir(x) = want [ fle—€) dw") (w > 0; » <n). 
0 


Hence, for —b < x < 2b, 
lyo?(x)| < C8 (5.5) 
by (3.1) and (3.2). Again, fora > 2b we have (x) = Oandf(x) = 0, 


so that 
x 


Iysr(a)| < [ \fl@w—€)||dw(E)| < CBe~* | e%|dwi(E)| < O8e~* 
xr—2b xz—2b (5.6) 


by (3.2) and Lemma 2 with o, = —c. By a suitable choice of C we 
have (5.6) true for all x > —b. 

Let us suppose that, for —b <x < some X >0, yx) exists 
and is unique for vy < m and (5.3) is true for v <n. If we write 
x = X+z2’ and apply Lemma 5, we see that y exists and is unique 
for x << X+C, andv <n. Hence, if y” ceases to exist (or to be 
unique) for any v < ” and any value of 2, (5.3) must have been 
violated for some smaller value of x and some v < n. For such pr, 
y”(x)e is continuous and so, if Theorem 1 is false, there must be an 
x > 0 and aA < n such that 


|y(ax)| = ee, |yE)| ee (—b CE< a; v <n), (5.7) 
whence |a(f)| << ee-E-) (0 << E<a; 8 < Ce). 
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By (5.4), we have 


Yew) = ya) I+ f Q(e—E,y) dw(E) 
and so, by (5.6), (5.7) and the properties of Q, 


eo < CBe-*+ | 0(ce-“2£) | du(E)| 
0 


< CBe-* +0 (eee [ e€ |dw%(E)|, 
a 


whence, by Lemma 2 with o, = —ce, 
e < 05+Co(e) < he 


for small enough e and 6 < 8)(e). This is impossible and so (5.3) is 
true for all vy < n and all x > 0. From this it follows at once that 


ly(a)| < ly) |+ j Q(x—E, y)||dwXE)| 


C(5+0 (€))e-* < ee-™, 


just as above. This proves the theorem for small enough « and its 
truth for larger ¢ is an immediate consequence. 


6. Extensions 

In the foregoing we have supposed that Q does not involve y™ or, 
more precisely, that Q satisfies a certain set of conditions expressed 
in terms of the first (n—1) derivatives of y. If we suppose 2 to 
involve y, we cannot readily apply our arguments as they stand. 
In Lemma 4 we should have to consider 


u(x) = | x(€,u) deo(e—€), 
0 


where w(x) = w(x) is of bounded variation but not continuous at 


x = 0. In fact, 


wnr(0+ )—w™ (0) = w™0-+) = k WHHE k,,(0) on 


The argument breaks down, as does the proof of (5.3), unless we 
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know that y™ is continuous. We can, however, apply our methods 
in an obvious way to prove 

THEOREM 2. If (i) 0, < —c < 0 im Lemma 2 and (ii) for allx > 0 
and for all functions y with continuous n-th derivatives satisfying 

lyY(a)|<C, («& >—b; v <n), 
Q(x, y) is a continuous, integrable function of x satisfying 
2(@,y)| = o() 
uniformly in x as B > 0, where 
B= Plz) — max |yé)|, 
<a—E<b 


0<z— 
ven 


then, for every positive « > 0, there is a 8 = &(e) > 0 such that, if the 
solution y of (1.5) satisfies a 5-set of initial conditions and if y™ is 
continuous for all x > —b, we have 

Iye)| < ee 
for all vy < n and all x > —b. 

Another type of extension to Theorem 1 deals with the equation 
(1.5) when Q contains linear terms with small variable coefficients. 
It is as follows: 

THEOREM 3. Theorem 1 is still true if the condition (5.1) is replaced by 

\Q(a,y)| < 0 (x)+-g(x)a(x), (6.1) 
where g(x) is bounded for x > 0 and 


fo) 


g(x) max J ef |\dw(é)| < (1—C)e” (6.2) 


von 0 


for some C > 0 and all x > some x. The last condition is certainly 
satisfied if g(x) > 0 as x >. 

This follows from two lemmas, whose proof I omit. 

Lemma 6. Theorem 3 is true if x) = 0. 

This may be proved by a fairly obvious variant of the argument 
of § 5. 

Lemma 7. If g(x) 1s bounded for 0 < x < xq and (6.1) is true, then 
for any «, > 0 there is a 8, > 0 such that 

lyM(x)|<g (vn; —-b<x<cnu) 


for any solution y of (1.5) satisfying a 5,-set of initial conditions. 


This is proved by the use of Lemma 3 and arguments similar to 
those of (17). 
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Summary 
THE problem dealt with in this investigation is the determination of 
boundary conditions in order to maximize a functional of both the 
boundary and the associated two-dimensional flow. The method of 
variation consists essentially in making a small local change at one 
part of the boundary, keeping the other conditions unchanged. The 
effect of such variations both in the physical and the hodograph 
planes is considered in detail for a number of problems, viz. solids 
of least resistance, solid of minimum mean hydrodynamic mass, and 
solids giving rise to greatest and least disturbances at large distances. 
The paper concludes with a discussion of a special type of problem 
where the solutions, which in general can gnly be determined for 
Laplace’s equation, hold also for the more general equation of elliptic 
type. 
1. Introduction: functionals of velocity in a field of flow 

A great amount of work has been published on the solution of flow 
problems starting from given boundary conditions and proceeding to 
the solution. For flows satisfying Laplace’s equation in any of its 
forms, these conditions are naturally given in terms of values of 
stream function or potential on certain surfaces. For compressible 
flow in two dimensions it is much easier to solve problems with 
conditions given in terms of the hodograph coordinates. In this 
paper, however, the problem considered is one which seems to be of 
equal importance, viz. the determination of boundary conditions in 
order that some desirable property of the flow may be as large as 
possible or an undesirable property reduced to a minimum. The 
problems which are in mind may be expressed briefly by saying that 
‘some functional of the boundary and the associated field of flow is 
to be made a maximum or minimum’. For example, to mention one 
problem, although a rather special one as will be seen later [§ 5], it 
might be demanded that an aerofoil in streaming flow should have 
the greatest thickness ratio for a given maximum velocity on the 
boundary, and so, of course, in the whole field. This problem has 


(Quart. Journ. of Math. (Oxford), Vol. 20, Sept. 1949, pp. 166-89] 
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been dealt with elsewhere but it is referred to for two reasons—it 
suggested the present work and it is typical in this respect: a func- 
tional of a geometrical character of the boundary curve, the thickness 
ratio, is being compared with a different type of functional, namely 
the maximum velocity in the associated streaming flow. If the 
boundary undergoes an infinitesimal variation, the thickness ratio 
changes at most infinitesimally but the velocity may be changed by 
a finite quantity. For example, if a small semicircular bulge is made 
on the boundary, the velocity is changed locally by a factor which 
varies from 0 to 2. If then the boundary is changed infinitesimally 
but without restriction on the shape of the change, it is not to be 
expected that a functional which involves the velocities will be 
stationary at a maximum or minimum. In fact, as will be seen later, 
it will be changed by an infinitesimal quantity of the same order as 
the variation, unless the latter is further restricted. The fundamental 
principle, in essence merely that used in all applications of the 
calculus of variations to maximum/minimum problems, will now be 
stated: If a functional of the geometry of a curve and the velocities (but 
not their space derivatives) of the associated hydrodynamic field is 
maximized by a certain curve, the functional is stationary for all varia- 
tions in which the physical coordinates AND the velocities are changed 
infinitesimally. It is not stationary for more general variations. In 
practice, a slight relaxation is made to permit finite changes in a 
region of a still lower order. No further attempt will be made at 
this point to justify this formulation of the basic idea of the calculus 
of variations. It is verified by the discovery that, for arbitrary 
infinitesimal changes in the hodograph boundary, only one equation 
results, whereas in the physical plane the variation depends on the 
shape of the infinitesimal change. 


2. The method of variation 

A method frequently used in applications of the calculus of varia- 
tions is to suppose a variation which is identically zero except in the 
vicinity of one point. This remark applies equally to the derivation 
of Euler’s equations and the use of variational methods in the trans- 
formations of the equations of dynamics and mathematical physics. 
An appropriate variational method for flow problems would seem 
therefore to be to apply a bulge to the boundary of such a shape that 
the velocity is changed infinitesimally, and to calculate to the required 
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order of accuracy the change of the functional due to the changed 
field of flow and the alteration of the boundary. In order to simplify 
the working and bring out the essential features of the problem the 
geometrical part of the functional is usually supposed to be just area 
(at one point merely a function of cartesian coordinates [see § 5]) and 
the main labour of these applications is in the calculation of the 
changes of the field of flow. For incompressible fields the technique 
applied is to map the boundary on a strip (first suggested to the author 
by Dr. W. H. J. Fuchs as the simplest procedure in such problems) 
and to take as variation a small ellipse in the plane of the strip. 
Exact expressions may then be given for the disturbed flow in terms 
of the original mapping function. Strictly speaking, the ellipse of 
infinitesimal ‘fineness’ does not satisfy the condition of producing 
only infinitesimal velocity changes, but these changes will be con- 
fined to regions which are of a smaller order than the ellipse itself.* 
The choice of the ellipse is purely one of convenience in reducing the 
labour of calculating the variations. An essential feature of two- 
dimensional Laplace fields is that such variations can be expressed 
explicitly in terms of the original flow functions. For three-dimen- 
sional fields and compressible flow this is not as yet possible.t For the 
latter and certain other problems, e.g. Stokes’s flow, a more general 
property of the variation may however be used to derive the solution 
in a certain special class of problems. These will be discussed briefly 
at the end of the paper, the main object for the present being to 
derive explicit solutions for two-dimensional Laplace fields. The 
corresponding step in classical problems is the formulation of the 
Euler equations for an extremal. It remains to state an elementary 
result which is useful at several points. 


Lemma. If the circle |\f| < p is mapped on a simply-connected area 
in the z-plane by z = f(C), where 


dz oe 2 bs 
dt =f (¢) in a 2 an (u/6) ? 


* Those portions of the bulge and the original contour on which the 
changes of velocity exceed any small quantity, given in advance, may be 
made as small as we please in comparison with the bulge by making the 
ellipse sufficiently ‘fine’. 

+ A generalization of the present method in which the whole boundary is 
varied in a suitable fashion avoids this difficulty and makes it possible to 
formulate equations for a wide class of linear-flow equations for which com- 
plete sets of solutions are known. 
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the a,, being complex numbers independent of yw, then the area enclosed 


in the z-plane is 
=. 2 
wofi-S} 
r 


The proof by integration of Z n round the circle is obvious. 


3. Bodies of least ‘dissipation’ and ‘drag’ 
It is convenient, and the reasons for the terms are apparent, to 


define d 
[to ds and | g(v) = ds 
y C 


as dissipation and drag, where C denotes the contour of an aerofoil, 

v is the velocity on the boundary, ds an element of the arc, x is 

measured in the direction of the stream, and the functions f(v) and 

g(v) are left arbitrary. The simplest problem of this type arises by 

making f(v) = kv and minimizing J/VA, where I = [ ke ds is the 
Cc 


linear dissipation and A is the area enclosed by the profile and 
circulation is absent. The square-root sign is inserted to make the 
ratio dimensionally correct, and the problem might be restated as 
that of finding the best distribution of a given A = Ag in order to 
minimize J. This problem, unlike most problems of this type, can 
be solved by an elementary argument as follows. 

Let the required aerofoil in the z-plane be mapped on the unit 


circle || = 1 by the relation z = f({) where 
f(t) = 14+ da,c. (1) 
The flow function is, with the usual symbols, 
w= $tip = UE+E). (2) 
7 . (so 
Then I= fe 3s ds = k{}, 


Cc 
and after insertion of the values at the limits for the upper and lower 
surfaces J = 8Uk and is independent of the coefficients in (1). It is 
this special feature that makes such a solution possible. In general 
I would involve the a,, in a complicated fashion. On the other hand, 


by the lemma, 4= the \a,, 2/(n—1)}. (3) 
The smallest value of J/VA is therefore derived by making a, = 0 
(n = 1, 2,...), and the solution is the circle. 
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Before proceeding to solve this simple problem by a method of 
variation it is convenient to introduce a certain notation which will 
be adhered to henceforth. Let z = 2+iy be the plane of an aerofoil 
and elsewhere referred to as the physical plane. Let € = €+% be 
the plane of the strip —2 < € < +2, € =d denote the position of 
the disturbance, and ¢, = €,+%, = pe'X the plane of the modified 
strip (see below). Further e is a parameter defining an ellipse; p is 
a small real quantity. Also ¢(B) is used to denote the value of ¢ at 
a point designated by ‘ B’ in each plane, and similarly for other points. 

The aerofoil being suitably mapped on the strip in the {-plane the 
complex potential of the undisturbed flow may be written 

w= Ul. (4) 
If the ¢,-plane is defined as follows 
C= d+l+pre*/ty (5) 
the small circle £, = pe‘X is mapped on a small ellipse enclosing D 
in the {-plane. 

The complex potential of a disturbed flow in the physical plane 
having the same boundary conditions, except for small elliptical 
bulges around the images of D in that plane, can now be given as 


w, = U(E,+p?/e,). (6) 

The varied J, which will be denoted by J,, is 
I, = 2k[$]g = 2kU [GIG [—1— w/a BE (@)]. (7) 
Now [<8 =4= (6,18 [1—pPe?/f,(B)o(@)], (8) 
where B, G denote the stagnation points in the physical plane, so that 
I, = 8UM1+p2(1—e2)/(4—d?)] + O(u3). (9) 
On the other hand, if Ay denotes the original area of the aerofoil and 


A, the area with bulges, 
|dz|? _ |dz|? 


A,= Ayt dnp dt Tar |e) +0(.4), 


de 
The condition that J/vA be stationary may be written 


8I/Iy = 48A/Ay. (10) 


Therefore 


(11) 


dz|” Si} 
2! D 


pe —ety(4—d) > 3 (1—etyutl fe 


di 


as e > 1, where suffixes 1, 2 refer to upper and lower faces of the strip. 
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That the circle is a solution is readily verified, and there is no other 
solution symmetrical about the axis of x. Thus putting 


C=Z4+23, Z=ee, 


we have ¢ = 2cos8, = 2cos8. 
Then (11) becomes, for e > 1, 
4 cosec?5 = sal La! f = "| eoseo%, (12) 
2A\|\dZ|,  |dZ}, 














and for the case of symmetry the only solution of (12) is derived by 
making |dz/dZ|? = 1. Also Ay = z, which gives the maximum value 
of A without any further calculation. 

A further point which is interesting is the consideration of (11) for 
the case when e is not nearly unity. The right-hand member is 
smaller than the left by the factor 4(1+e?) and this proves that 
5I/Ig > 48A/Ay, which indicates that the solution makes I/VA, a 
minimum. The fact that this necessary condition for a maximum 
or minimum can be demonstrated without the necessity of calculating 
the coefficients of higher powers of yp is of great utility because the 
computation of even a first variation is in more difficult problems 
sufficiently tedious. That the variation is not stationary for bulges 
with 0 < e<_1 is in accordance with the principle stated in the 
introduction. 

The nature of the method having been indicated by the simple 
example, it remains to carry out the analogous process to find an 
equation corresponding to (12) for the more difficult problem of 
minimizing R®/VA where 


R® — | m(Z) ie: 
Cc 


In this problem there is no such simple expression as (7) for R, since 
the latter involves a complicated function of the A,,. In view of some 
results to be discussed later [§ 5] it is useful to note that the general 
character of the solution of this problem is that of a fairly fine oval 
shape. In a general way this may be seen by comparison with the 
circle which solves the preceding problem. Both the factors v? and 
dx/ds are too large over the middle portion of the circle, and so the 
minimizing curve is generally flatter than the circle. Set 


dw|? 


(2) 
R lz dx, 
o 
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and to save various complications a solution will be sought which 
has symmetry about the axis of « drawn parallel to the stream and 
making the ratio R®//A stationary. Here w is given for the disturbed 
flow as in equations (5) and (6), the relation between the z- and 
¢-planes containing the unknown functional dependence. Now 


len. 12 den ft 2 
er lz = U*|(1—p?/23)|2/|1—e%2/E3|%. (13) 


at | ~ \dt,| /\ae, 
io alt Do(e) } i 


From equation " 
where the C, are binomial coefficients and positive. 
For all ¢ on the strip and outside the ellipse around D, 




















|2eu/(f—d)| < Qa- =b6, say. 
Then (14) gives |we/f,| < [1—(1—b?)#]/6 = e. But 
po | | ae (1—e? )ye 9.9 2) 1—e? 
|< 14-4 2e2/C2|\"] <1 ~ 
aah s d+ fe [ + > |pPe*/2i| ]< 7-3 


It readily follows that, for all e < 1, the result of substituting 
(14) in (13) expanded as a series is absolutely and uniformly con- 
vergent with respect to { in the range of integration. The problem 
is now reduced to the evaluation of 


my = [ F(X (E—d) dé 


along the part of the real axis indented at D. Let 





A= + ¢), 
Then said 
Pace §) d& pe’ 
als eae a  (€—d)2n - 
— a(§,d 2N 1 —d)-2N+ 
oe dex ) dé+F(d)u |-aai¢ d) vel 4 


+P (dyer |- — (gaye . (15) 


For N > 1 there is no term in p?, and 
Ke’; 


ty = 2F(d)p24/(2N— 1)A2N-1 3 2F(d) | cope ——_. dé+ O(p3 ), 


pen 
(€—d)? 
A+d 
where A,(§,d) = F(é)—F(d)—(§—d) F’(§—d). 
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Let R, and R, denote the values of f v* dx taken along the real axis 
é 


and the elliptical indentation respectively and write f(é) for the 


boundary value of da 
“(as) 
Then 
2 d-~A 
k=[ f+ f |pon—2a—emie—arjt2r@ [Sa 6) 
d+X -2 d+A 





y2NOy _ [ 1—p2/C2 \? 2y?(1—e2) 
where >= DE a = (ie — ae . 
The second term can now be evaluated as a definite integral where, 


for the first part, it is convenient to change the integration to the 
(,-plane. The result is 


1 2 1—e? 1 
afta) [ Sag =: a f+ Sloe (*4)]. 


d+A 





Again | f+ re pas fred f(é) dé—2f (d)-+ 008) 


d+aA -—2 


o-oo fT 
stl ood | ieee at 
The final result is 


__¢2)2 l 
n= [16 £) dé—2uf(d \|e- — So" toe(**}]+ 
—2 


3—d d) 
+2p%(1— of £8 —f' @los(5=4)— Ses | (17) 


The existence of all derivatives of f() is assured by the fact that f(£) 
is a rational function of the derivatives of a function of a complex 
variable in a region free from singularities, the ends of the strip being 
for the present excluded.* 


and 








1—e? 
= = 4S S)r@+ 




















* The difficulty at the ends is only apparent and does not arise if the 
strip is mapped on a circle. 
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The evaluation of R, to order »? must now be carried out, and it 
will be found, as might be expected, that the terms in p cancel. 
Indeed the appearance of a term in is merely a consequence of the 
rather artificial method chosen in the calculation. In the example 
of the next sections this is avoided by a special device. 
dw 21df,|? 
Ta = 


2 








Now R= | 














dt, 


and, the length of the path being of order y, the integrand must also 
be found to order y». To this order the following holds 


Ax+iAy = (F ---$ %) AC 


-(F at’ aE) —(-+e)sing +i(1—e4)o0s x}u dx 


for changes along the bulge. 


- at (a d (dt 
_— az ~ (a), [lee 2 
dz\-* = (dz\- d*z 
(a = (zi), ‘ * wag i). 
where AC = p(1+e?)cos y+ tc(1 —e?)sin x. 


It is convenient to put 














Then 
de he = de i {1—2p[(1—e?)a cos y—(1—e?)B sin x ]}+ O(u?) 
\d¢| dl\4 
: dx . dy ; : 
Again — S41 = (we bingda + (rge tinge AL+ OU). 


dé" dé 
From these equations it appears that 


F 4sin*y dx 
Ki ye | I ie 
. | ; 1 —e?)? cosy + (1-+-e?)* sin? y ‘ 
0 
x {1—2yu[a(1-+e?)cos x—B(1—e?)sin x ]} x 
x {ae(1-+-e?)sin x-+-yg(1—e?)cos y+ 
+-age(1-+e4) sin 2x-+-yee(1—e4),. 008 2x34]. (18) 
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r sin®y dx ae 
D 2(1-+e?) 
D = (1—e?)? cos*y+-(1+e?)*sin?y, 


- sinty dx = te apes +5] 





D 1+e?)’ 


F sin®y cos 2xdx__—_—a(1—e?) 
D ~  &(1+e8)” 


0 
Therefore 


] 1—e?)? 1 
R, = “ae 1— OS hog(**)] + 


al 


+amp(1 [(1—e*)yeet-aye(1 +e?) +-Bare(3—e*)]g. (19) 


As in the preceding example the — a area is 
5A =p tn(1—et)| (20) 


The equation which connects the PE sg values is then found by 
putting the coefficient of ~?(1—e?) in R,+R, proportional to that 
in 6A and letting e > 1. We obtain 
1 dt|? dz|? 
-7 = |— paaed 21 
aS(6) =|72) Weel #7 (21) 
where the linear operation 7) is that of equation (17).* 

It is apparent that a similar result holds for any J or R and that, 
moreover, an equation can be found to determine the form of f(v), g(v) 
in order that a given curve, e.g. the ellipse, may give maximum 
values to I/VA, R/WVA. 


4. A problem due to G. Pélyat 

The problem of this section has been solved by G. Pélya [2]. It 
may be stated as follows: T'o determine the shape of the aerofoil which 
in unlimited streaming flow has the least mean hydrodynamic mass 
(with respect to all directions of the stream). That the solution, as 

* In the plane of the circle the operation is equivalent to taking the 
derivative of the conjugate Fourier series. In terms of a more general 
variation of the boundary this may be given a simple physical interpretation. 


+ The result of the section follows very readily, as shown by Pélya, from 
the Lemma of § 2. 
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proved by Polya, is a circle is, of course, very plausible on grounds 
of symmetry. The method of the preceding section will now be 
applied to derive a functional relation which is satisfied by the 


circular boundary. 
With the notation of the previous section let the plane of the 
modified strip with circular bulge be mapped on another strip in 


the t-plane by - 
, , t= 0+ptr?. (1) 
The flow function for a.stream inclined at « with the positive 
x-direction and with the slightly varied boundary conditions is 
w, = Utcosa—iU sin «/{(t—t{ B))(t—t(@))}. (2) 


In terms of ¢,, 


l 
w, = Ucosa|tp+¢,—f,(B)+ a: ~~) |- 
—iU sin ay/{[,—2(B) [21 —41(@) 1 —?/21 Sy (B) JL —p?/L, 23 (@) J}. 
The ae mass M is Pie het 27’ = Mu? where 
2T = p | (;2+¢;2) dS 
and ¢’ is the potential for the disturbance produced by the moving 
cylinder. This may be expressed as a line integral as follows, where 
for convenience it is supposed that p = 1, 
or = [ ods. 
J ve (3) 
The integration is along the boundary, and clearly the plane of 
integration is immaterial. Let M denote the mean value of M for 
all « between 0 and 27. Then 
27 — MU? 


2a 


U2 


2a 


| [(A—«x)cos a+ (B—y)sin a} E cos w+ a sin o| tsi 


0 C 


= 10? | (4-2) F4-(B-y cubs, (4) 


on 


é 
where A is the real part of 


1 
tpt+O—o,(B)+p? (-- nm) 
B is the real part of 
E(B) —&) (G.—Sa(@) (1? /bs S(B) (1/0 @))}- 
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Briefly this may be expressed as 
2M = i F(t) ds; = | (Fy+5F) ds; 
C é 


= | RO d+ [ aF dst, (5) 
Cc Cc 


where F' is the function in (4), K is the corresponding function 
for the undisturbed flow, and the accent denotes the varied contour, 
viz. the strip with a small elliptical bulge. Now the function of { 
from which F, is obtained has no singularity at the bulge, and so 
the first term of (5) may be replaced by the integral of the original 
function over the undisturbed strip together with a surface integral 
over the bulge. The latter may be evaluated by Green’s theorem. 
It is convenient to consider the portion above the axis of £. Thus 
let A and a be regarded as point functions of _ n) and 


| (Vxr){V(Ay—ax)} dS; = | a —2) = = dy (6) 


Now Va = (x,%,), 
V(Ay—2z) = (Age—2z, Aon—%y), 

and in evaluating the surface integral to order ,? it is sufficient to 
calculate the derivatives at d. 
Now (Va){V(Ay—2x)} = (Ag—axg)ug+(A,—2,)@y, (7) 

(Vy){V(Bo—y)} = (Be—yel¥et+(By—Yn)Yq (8) 
The contribution of (6) to order p? is found by multiplying the area 
of the semi-ellipse by the sum of (7) and (8) calculated at D for the 
undisturbed functions which are 


A= é, B= W(4—&) ; 
we got Anat —e)[ yy tts — 2108 +98) (9) 


The term of order 1 in (5) is the undisturbed value of 2M = 2M, 


where ef [ (A— 2) + (B— y) alae. (10) 
—2 
The second term of (5) must be evaluated in two parts, the first 
along the real axis between 2 and —2 without the portion 
d—A<&<d+a 


3695.20 N 
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and the second along the upper half of the elliptical bulge. By a 
calculation analogous to that in the preceding section it is found 
that the first of these is 


(¢)—g(d) d 
eye dl+g(djlog | (11) 
Jeg acsnaind = 
where g(é) = ox d(4—£€) oy 


on 4—d? én’ 
and the error is O(u3). The last statement follows from the double- 
series expansion for the varied potentials in terms of { [cf. equation 
(14) of the preceding section]. 
The second part along the bulge requires a knowledge of the 
derivatives to order pw and, to this order, 


8A = R(1—e*)p?2/L,, 


d 1 
= —Rul—e*),(4—2) "=. 
5B Rp*( e*),/( T-@L 
Also 
de ded pre ; . 
rnd i (—ixy-+y,)(1—eu2/£2) = (am,-tiyme-®. (12) 


The required integral is 


[er ds = | an? 2H) as 
-nfo-onfl2)-(Bhctal 


which after using (12) to find the derivatives in terms of x becomes 


oe dy da d 
=a-emttel)+(Gha—wle 


The contribution of the upper part of the strip is then found by 
combining (9), (10), (11), (13) to give 
2M, = 2M,+852M 


2 
== | a Ad dé + p?(1—e*)[T 9(€)]+ 
-2 


+p2m(l—e (a2 +h) —H* m(1—et)(x2-+y2)q, (14) 


e— dF) 


where 7), is the linear operation indicated in (11). 
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On the other hand, to order py”, the change of area is, for the upper 
bulge, 8A = $mp*(1—e!)(22+y2). (15) 


It is now easy to verify that a symmetrical solution is given by the 
circle, for on putting 


f= Z+2Z-1 2= 2 2(d) = e® 
and, evaluating the right-hand member of (14), we have 





a 2(3—e%)(1—e2 
nighe [i+ a | si 
H4o+S4) = do[l+eA1—e\(1+e4/45in%8}, 27) 
5M 5A — = 
a or 2(1—e?)[ (1—e?)/2 sin®8], 


and so, for e < 1, the ratio M/A is increased and, for e = 1, it is 
stationary. The circle is therefore a solution and it gives a minimum 
value for M/A. A rather better discussion of the solution can be 
made by transforming (14), (15), and the corresponding results for 
the lower side of the strip to the plane of the circle. Let the (unknown) 
aerofoil in the z-plane be mapped on the circle Z = e and suppose 
that 








dz : 
eet 18 
& = H(0)+iG(6) (18) 
on the circle. Then the boundary values of x,,y, are given by 
: H+iG)e® 
ty tiyg = EE (19) 
1 , , 

g(é) = Fan®s [ H(@)(sin®d cot 9—cos 8 sin @)— 
— G(@)(sin?5-+-cos 8 cos 6)], (20) 

T,g9(é) be} fom (0)+-H(@)cos @ cos 8} as| 

tal | sino sing SO—8O) d+ cos j G(6)sin 6 as| 
ya — cos 


0 





G(8) 1—cos8 
——. log ° 
2 sind 1+cosd 
+t The order of the variation is u*(1—e*). It is stationary because of the 


second factor 1—e®. In other words the terms in (16), (17) cancel very nearly 
ase —> 1. 


(21) 
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For a symmetrical solution H and G are conjugate Fourier series in 
cosine and sine respectively, and in the case where symmetry is 
lacking the sum and difference of the values on the upper and lower 
surfaces again give rise to conjugate series. 


After a suitable reduction, the variation of M, is found as 


p? SF [aha H +e) H+ G+ H+ @B))* 
The variation of area is 


4 u2(1—e?)n[ H?-+ G?+ H3+ G3]/4 sins, 


where 2H = H,-+-H, and the suffixes 1 and 2 refer to the upper and 
lower surface. It appears therefore that it is not very easy to prove 
that the solution is symmetrical by this method but, if symmetry 
about some line is assumed, the equation of variation demands 


H?(d)+-G?(d) = 2k? 
for all points, where k is some constant. ad the unique function 
regular in Z > 1 which takes the value | f’(Z)| = 1 on the circle is 
1. Therefore H+iG is constant in whole dig |Z| > 1, and the 
solution must be the circle. 


5. Variations in the hodograph plane 

When for the purposes of a problem the upper limit of the velocity 
is specified, it is often convenient to work in terms of the hodograph 
plane defined by Q = i dw/dz. The factor i is introduced merely for 
convenience as mapping the upper half of the physical plane on a 
region where amp Q lies between 0 and z. For a given flow function 
in terms of Q the coordinates of the physical plane are given as 
functionals of the boundary in the hodograph plane according to 


the relation 
tliat a 


i 
A local variation applied to the hodograph boundary changes the 
whole physical boundary of the flow and any infinitesimal bulge 
corresponds to infinitesimal changes of velocity in the physical plane. 
According to the general principle of the introduction, if a functional 
of the flow involving the velocities, but not their derivatives, is 
maximized by a certain hodograph boundary, this functional must 


(1) 


* Where h, is the first coefficient in the cosine series H. 
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be stationary for all changes. It is therefore very desirable to verify 
analytically for these special changes in the hodograph plane that 
the variations in the problem considered are independent of the 
shape of the bulge. The proof for the coordinates (x,y) themselves 
depends on the fact that z is a holomorphic function of Q. The 
more interesting problem is to show that the area enclosed in the 
physical plane has the same property. This is proved by the device 
of splitting the field into two parts, one including the bulge and 
extending to a distance of order vy, and the other the remainder 
of the field. 

if w = w(Q) takes real values on a certain curve in the hodograph 
plane and has the correct behaviour at the singularity iU, where U 
is the stream velocity, the complex potential of another flow having 
the same hodograph boundary except for an indentation at D is 


w, = w+e/[w—w(D)]. (2) 
The varied physical coordinates are given by 
a= [Fae /fo—wiD)p (3) 
é 


In the neighbourhood of the singularity iU, w becomes infinite, but 
by appealing to the physical plane it is easy to see that the integration 
for the variation occurring in (3) may pass through this point. It is 
therefore always possible to calculate a change of « and y between 
two specified points in the hodograph plane by a path along the 
boundary but avoiding D. It is then evident that, except in the 
vicinity of D, (2) represents the complex potential due to any bulge 
to order «?, where « is the order of the linear dimensions of the bulge. 

For example, let x be the chord in the direction of flow for a 
symmetrical shape having a simple representation in the hodograph 
plane and y the maximum ordinate. Also let M be the point where 
the hodograph boundary cuts the imaginary axis. It is convenient 
to map the flow on an auxiliary strip and it is found that 


M 
e { cosO(é) dé 
w J q(é)(é—d)” 


‘ff flees a)" 


Ay = (4) 
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where Q=ge®, f=é-+in. 
Except for a constant factor, these variations hold for all bulges. 
It is also noteworthy that, for such simple curves, sin@ > 0 and, 
if d lies in the first quadrant of the hodograph plane, cos 0(€) corre- 
sponds to the second quadrant and is negative. The signs of Ay and 
Ax are in fact essentially negative and positive respectively. If, 
therefore, it is required to make F(x, y) stationary and so 
F,Ax+F, Ay = 0 (6) 

and if the partial derivatives of the function F have opposite signs, 
there can be no solution. It was by a similar reasoning that the 
author, by taking F = y/x first discovered that no extremals in the 
hodograph plane exist for the problem of making thickness ratio 
a maximum with a given upper velocity. The correct answer to this 
problem, as had been anticipated on other grounds, is derived by 
letting the hodograph boundary take extreme values [1]. The result 
is the same for the whole class of problems where F, and F, are of 
opposite sign. It also appears from a consideration of (4) and (5) 
that, when F,, and F, are of the same sign, e.g. / = ay, the solution 
has singularities at d = 0, d = 2, corresponding to the ends and 
the middle of the aerofoil. It is therefore probable that no complete 
solution exists for any such problem, but that further restrictions 
enter and the extremal condition applies over a limited arc only. 

Before entering upon the calculation of the variation of the func- 
tional A = area = [ ydx 

é 

it is convenient to state a problem which is associated with the 
vanishing of this variation: it is required to maximize the ratio D/A 
where D is the disturbance of a body at great distances and the velocity 
is restricted on the boundary. It is readily seen that the restriction 
is necessary if a definite result is to be obtained. Thus in the notation 


A = m{l— > a2 /(n—1)}. 
The complex potential for streaming flow is 


w = U(t+p2/2), 
df = C—a,p?/6-+0(E-). 


of the lemma 


dz 
dg 


where z= [ 


At great distances 


w—UL = Up(1+a,)f4+O(f-) = Up2(1-+a,)2z3+ O(e-). 
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Therefore D = yu?(1—a,) is a measure of the disturbance and it is 
also necessary that 1+-a, > 0. For a proof that this is always true 
a slight rearrangement of § 2 [1] is sufficient. 

Now A/D = x{1—a, — > a2/(n—1)(1+4)]. (7) 
Then A/D has a maximum 27 for the limiting case of 

a, = —1+0, a,=0 (n> 1), 
a thin ellipse along the direction of flow, and it has a minimum of 
zero for the case a, = 1, a, = 0 (n> 1). This second result is, 
however, trivial because the same is true for any very thin body 
placed across the stream, and the velocity increases without limit 
as D/A — ©. 

The problem with the restriction on velocity is solved in the 
hodograph plane by equating the variation AA to zero since, for 
changes of the hodograph boundary, the conditions at the singularity, 
corresponding to the physical flow at great distances, are not affected. 
It is evident that there is a certain portion of the hodograph boundary 
lying on the circle |Q| = M, where M is the given limiting velocity, 
and the equation AA = 0 determines the extremals only. 

In order to calculate the variation of area it is convenient to 
introduce the notation of preceding sections. Let the bulge in the 


plane of the strip be derived by transformation of the circle {, = pe*x 
in the modified strip. Thus 


[= d+ G+ wel th > (u/¢,)"—", (8) 
3 


where the a,, are independent of yu and the series converges uniformly 
in the region |f,| >. The complex potential of the varied flow is 


w, = U(%,+p?/%,). (9) 
In the region |f,| > vu the following equations hold: 
C—d = C+ p%e?/f,+ O(n") (10) 
dw _ U—w/) 
a ~ iO — Ya, wr A 
= U[1—p?(1—e*)/G3+ O(u/f,)*).- 
C, = (—d+pO(u/¢), 
1 1 


at apa} 
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dw re p*(1—e*) ( a y (+)] 
sil >) PlF}I- 
Therefore dt | ({—d)? . t—d a 


= Job cau Hea) 4 
© = [oer tea) Ce) 
The integral of the second = m is of order 
Op'/(C—d)® = O(u!) = o(p?). 


Therefore, to order p?, 


From (1) 





/(1—e?*) 


where z, is the value for the ae Phi If the path of integra- 
tion is started at € = —2, then for points lying between —2 and d 


g 
—%= ~ing1t.—e] | pitt 
—2 


] 
— Pa 


See d— 
t+ycq}+2 (log \3 } 


U af . ' 
where P = —, and — is continuous at {=d. For points 


Q (¢—d)? 
between d and +2, 


—iz = —iz—p?(1—e?) x 


é 
«Leo halp at tal + rene 
; (14) 
For points on the semicircle |f,| = u' similar expressions hold and 
in this case the accuracy required is only ,! since the path is of length 
O(u*). Therefore, to this accuracy, 
y =yotReO a. 

Also f—d = pte*x[1+O(p)]. 
Therefore Y = Yotpi(1—e?){C(d)cos y+ S(d)sin x}, (16) 
where P= C-+48. 

Let A’ denote the area enclosed in the physical plane of the varied 
flow by the transform of the strip in the ¢,-plane indented by 
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\f,| = pw. This is not the physical boundary for the disturbed flow. 
Also let Aj, A} be the contributions for the straight portion and A, 
that for the curved. From (13) 


dx p?(1—e?) 
— = §,| 1——____-]. 
dé | (€—d)? | 


Therefore 


ay eenere'= | 

x {Sy}{1— wey (E 9 a, 

= { voSae—pr1—ery [ Wo J w(1—e)Sy(€) x 
1 1 


‘ 
A, C 1 
2 d&t—C@)|— d 
«| | Tape gaat ga} + OAs =il| 4 
—2 
Similarly 


at a §,dé—p? vol) Sol fea 
Ay J mesa pl 0) fe [ €2)84(£) x 


é 
As es. cal l . , 
«| | eae @\~ gtaalto" (d) og|S— ‘|+S@] dé. 


2 
4,+4, = oe | water ae + 
1+2 =.“ 
1 


nt 2 


pee ‘|| 


—p*(1—e2)78'(d @ f S\(é) dé, (17) 
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A, having the same meaning as in § 3, and A, denoting the difference 
from the value of the function at D. Again, on the curve |f,| = p?, 


(da, .dy\— w_-o\fy_ w(l—e) 
—i(E+iF) = is/ Ta | (18) 


Therefore 


-_ r{ —iS)) IF a toy | x 
x [Yo tv! (1—e?){C(d)cos x+ S(d)sin x}| df 
2 
sas | (Co—So)¥o at —p(1—e*)[(Co4—i8yye)e — 5) + 


+ (Cy Yo—tS Yo)a 7+ 7(1—e?) | (Cyo—t8p) x 


x (C, cos y+S, sin x)e% idy+O(u#), (19) 
where I denotes ‘the imaginary part of’. The suffix 0 refers to the 
value at that point for undisturbed flow. Combining (17) and (19), 
the terms in p! cancel and 
A’ = At Ay +A; 

i 2 
= Yo dxy—p*(1—e?*) | 73.¥05)) + C(d)T(So)— | SoQ)E7s C(f)| d¢— 


1+2+8 -2 
. 
=8'(@) | 88) d&—Cyyon+4(8+89)| 
d 
= [ Yo dty+AA’. (20) 
1+2+3 
The 7'’s are used as abbreviations for the linear operations involved 
in (17) and (19), and the result so far is just the same as if a simple 
doublet had been used having strength »?(1—e?). It remains to 
show that the variation of A is actually the coefficient of u?(1—e?*) 
in (20). First of all it is to be noted that the expression in (20) exceeds 
the derived variation by the area corresponding to the annular region 
p< |¢,| < pw? with the disturbed flow functions. On the other hand 
the first term of the right-hand member of (20) exceeds A, by an 
area corresponding, with the undisturbed flow functions, to the bulge 
together with the annular region. It remains to calculate these 
changes using the double integral for area 


[fon J SEBan= tes = [el ate 
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where dSg is an element of area in the hodograph plane. It is also 
evident that an expression of the form (21) may be evaluated in any 
plane. Then 

al 


Ai +A,+A) = Ap+AA+ alg dS, 


where w, is given by (9) and the region of shite is the annulus. 


Again 
dw 
f ey Say Jala dS;+ Sela 2 wi" a8, 


1+2+ 

the first pine (B) being over the bulge in the ¢-plane and the 
second (A) over the annulus in the {,-plane. Then equation (20) 
becomes 


1 |du 
A,+A og taseail 
ot i Plat 


1 dw|? ’ 


To order p?2, since there is no singularity in re undisturbed flow, 


dw |? ., ‘ a2 |U? 
jatze-~|"-- pile 
B 


by the lemma, and 


~ @ dw dw|? 
, ace 2 o4 n machen’ | 
AA—A’A = Th (1 _— > a i)- fe | al, ra | as, 
(22) 


2 
From (9), aioe 1 (2 wait F) : 
1 1 


ron Ene ofi-He- Ball] 
3 


Let ; dS;, = pdpdx. 


Then 


1 (|dw,|? |dw | 
—{|—— dS 
iF | a, - lat . 
U? 2p* cos 2x 7908 2x pt -| ‘m ay, —n)i anu | 
pes as. 1 pmtnein— ix S 
=| fie ==[1— p* TA a poe = mane” or y p" 
(23) 
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The term 


2 
[female peat 


q 
u* (m—n)i ‘ad coca O 4+m+n —m—n+1 
“thao ome pdpdxy+O | pitmtnp p 
bw 
= 0+ O(umintip-m-—mtt),_, = Oly). (24) 
The term 


2 


pw 
wo (HN dp = 2an™ | pad = L 
q°\p ¢ : gq n—1 
Be 


from which it follows that to this order (23) is* 


2 
mt 1—e— 2 me i} 
3 


Finally the right-hand member of (22) vanishes and 

AA = A’A. (25) 
This completes the verification that only one equation results for all 
possible shapes of bulges. 

The equation derived by making A’A vanish in (20) determines the 
extremals for the problem of maximizing D/A subject to a given 
velocity restriction. The solution is thus in principle reduced to a 
problem of the ordinary calculus. 


6. Compressible flow problems 

In the above work the reduction of various maximum and mini- 
mum problems has depended on the fact that the disturbed flow 
functions could be represented explicitly in terms of the original flow. 
This cannot, as far as the author is aware, be achieved for com- 
pressible flow or fields satisfying differential equations other than 
Laplace’s. There is, however, a rather special class of problems 
which may, broadly speaking, be characterized as those for which 
the maximum or minimum solution is derived by taking the extreme 
boundary of the permitted hodograph region. In the previous section 
it was observed that Az and Ay were for a certain variation always 
of opposite signs. A closer inspection reveals that this depends on 
two features: (i) the continuity of the streamlines due to the bulge 


* The error in taking g as constant is of a higher order than p?. 
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disturbance, (ii) the condition that z was a holomorphic function 
of Q. The first condition is equally true for subsonic flow whilst the 
second is merely a way of saying that dx and dy are perfect differen- 
tials of the hodograph variables. This, however, is the content of 
the equations of compressible flow. It therefore follows that if F,, F, 
are of opposite sign the extreme hodograph conditions permitted 
will be attained. The simplest condition is to confine the solution to 
a semicircular region giving in the case of F = y/x the aerofoil 
problem of reference (1) a flow due to Riabouchinsky. The result 
depends, however, on the assumed limiting conditions and is not 
necessarily that solution. 

To conclude, a very simple problem concerning circulatory flow 
can be solved by a similar argument. It is required to dispose of a 
certain area so that, in purely circulatory flow of a given vortex 
strength, the velocity on the boundary is as small as possible. If this 
velocity is not constant, let a small bulge be made at a point of least 
velocity and [cf. 1] the area is increased and the velocity decreased. 
Then reduction of the whole figure by a scale change to its original 
area does not change the velocities. Therefore the velocity is constant, 
and the solution is a circle for incompressible or compressible flow or 
any field satisfying an equation of elliptic type and with the same 
boundary conditions. 


REFERENCES 


1. A. R. Manwell, ‘Aerofoils of maximum thickness ratio for a given 
maximum pressure coefficient’, Quart. J. Mech. and App. Math. 1 (1948), 
365. 

2. G. Pélya, Proc. National Acad. Sci. 33 (1947), 218-21. 

3. O. Bolza, Lectures on the Calculus of Variations (New York, 1946; reprint 
of 1904 edition). 











A NOTE ON S-FUNCTIONS 
By H. O. FOULKES (Swansea) 
[Received 8 November 1948] 


1. THE characteristics of the classes of the symmetric group of order 
n! in which all the cycles are divisible by r can be obtained from the 
character table of the symmetric group of order (n/r)! by expressing 


S-functions of ont 
v1; €X1,; € Ly se005 € v7; 


Sig, Chay OC Mes.<-, Chg, 


(1) 


2. r—1, 
vy, EX; , € Lgseeey € x; 
in terms of S-functions of 


Wy, h,.00 Xe, (2) 
where « is a primitive root of 2” = 1., In this note I point out a very 
simple way of expressing S-functions of the first type in terms of 


those of the second type. 


2. Any S-function of a set of numbers a, ag,..., &» is a polynomial 
iN 81, S9,..., Where s, = }a?. Consider those S-functions {«} of (1), 
corresponding to integral values of x, including the case {0} = 1. 
Unless p is a multiple of r, s, will vanish. If S, is the product sfs},... 
associated with the class p, and h, is the order of the class, then 


1 
fg => > 4,5, 
Pp 


and so {x} will vanish unless there is at least one class in which each 
of 1, m,... is divisible by r. Hence, for {x} to be non-zero, « must be 
divisible by r. Corresponding to each class r*(2r)”... of the symmetric 
group of order «! there will be a class 1*2”... of the symmetric group 
of order («/r)!. If p, and p, denote these two classes respectively, thent 


x! («/r)! 
h, = a * = —— 7 
" pty! (2r)vy!... 17x! Qvy!... 
! 
and so i a oe Fn. 
Pr («/r)! pEt+Y to 
+ (1), (3), 143-6. t (3), 40. 


[Quart. Journ. of Math. (Oxford), Vol. 20, Sept. 1949, pp. 190-2] 
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Also, if z, denotes the sum of the pth powers of (2), then s;,, = rz;, 


and 
{x} = DD h,, #8... a — ps) h,, 2%... 
= fairy, 


where the superscript denotes an S-function of (2). 


3. Any S-function {A,, Ag,..., A,} can be expressed as a determinantt 
in which the term in the pth row and qth column is {A,—p+4q}. 
When {A,, Ag,.--, A,} is an S-function of (1), the determinant will have 
zero elements whenever ’,—p-+q is not divisible by r. The other 
elements will be of the type {x/r}®. Hence to express any S-function 
of (1) in terms of S-functions of (2) all that is necessary is to write 
down the determinant |{A,—p+q}|, replace every element of the 
form {rs} by {s}”, and replace all other elements by zeros. For 


example, {5421} = | {5} {6} {7} {8} 
3s 6 6 © 
Oo 8 8 8 
| {0} {1} 

When r = 3, 

{5421} — | {2}) : 

| ae , : {2}) 

| 1 ; ; oor 

| 1 


~ (31944219, 
Similarly, when r = 2, {5421} = —{37}, and 


when r = 6, {5421} = —{2}©—{1}, 


4, The same results may be obtained by expressing {A,, Ag,..., A,} 
aS |@,,-pig|> Where {y,, Me Ht is the partition conjugate to 
A,} and a, = {1*}. Every element of the form {1’%} is 
replaced by {1°}, a negative sign being attached when r is even 
and s is odd. The other elements of the determinant are replaced 
by zeros. 
t+ (3), 88. 
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5. S-functions of (1) and an additional number ¢ can be dealt 
with in the same way.} If fA, Ag,..., A,; 2}, or {A; 2}, denotes such an 
S-function, then it is known that 

(A; Y= =a {A} +[ DI IE+[ DS Popalet]o?+-. 
where {v}{u} = > 9,,{A}- The coefficients of the various powers of ¢ 
may each be reduced by the method of § 3, or {A; ¢} may be put into 
the form |{A,—p-+q; ¢}| and each coefficient in the polynomials 
fie; CO} = {«}+ fe —1} 0+ {«e—2}02+...+¢* 

reduced. Thus {321*; £} becomes, for r = 2, 
| \{1} — {2} 4+ £1 }2Z24 74 {2}2¢ +-{1}° 2734 75 {3} + {2} O72 4 (O44 76 


‘.. {1}¢ 142 {1l}Cf+28 {2}@) 4 £]}@724 64 
1 ¢ (1y+¢8 
a | | 


= —{19}@f—{ 12}(2)73, 


I wish to thank Professor D. E. Littlewood for some comments on 
this note. 


+ The results given in (2), 197-8 and (3), 135-6 are not quite correct. The 
indices of { need to be adjusted. 
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